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Abstract 

The  Aeronautical  Systems  Division  Infra-Red  Emission 
Prediction  Program  (ASDIR)  is  a  large  computer  program 
compiled  in  1974-75  using  "state  of  the  art"  routines  to 
predict  hot  surface  and  plume  emission  from  a  geometrically 
modeled  target.  The  atmospheric  transmission  segment  of 
the  program  included  only  the  major  gaseous  absorbers  (^0, 
CC>2 ,  and  with  no  aerosol  attenuation.  To  improve  the 
flexibility  of  the  program  in  making  more  accurate  predic¬ 
tions  at  low  altitudes,  an  aerosol  transmission  computa¬ 
tional  scheme  was  written  for  the  ASDIR  code. 

Using  the  basic  methodology  in  the  Air  Force 
Geophysics  Laboratory  (AFGL)  transmission  model  LOWTRAN  4, 
a  subroutine  was  written  for  direct  inclusion  into  ASDIR. 
Curve  fits  of  the  following  LOWTRAN  parameters  were 
performed  which  simplified  the  aerosol  procedures  and  mini¬ 
mized  the  additional  computer  storage  and  time  requirements 
in  ASDIR:  (1)  assumed  aerosol  particle  densities,  (2)  equiv¬ 
alent  sea  level  path  lengths  for  vertical  paths  -  AV(z) , 

(3)  total  sea  level  extinction  coefficients  for  an  "average 
continental"  and  "maritime"  aerosol  models.  Results  of 
the  new  subroutine  (called  HAZE)  compared  very  well  with 
LOWTRAN ’s  aerosol  predictions.  A  logic  change  was  included 
in  HAZE  to  automatically  select  a  maritime  aerosol  model 
under  low  visibility  conditions.  This  programming  change 
was  prompted  by  results  of  a  recent  Grafenwohr,  FRG  trans¬ 
mission  study. 
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ADDITION  OF  AN  AEROSOL  TRANSMISSION  MODEL 
TO  THE  AERONAUTICAL  SYSTEMS  DIVISION 
INFRA-RED  EMISSION  PREDICTION  PROGRAM 

(ASDIR) 

I .  Introduction 

In  recent  years,  applications  of  infrared  technology 
have  multiplied  in  the  Air  Force  and  other  services.  Of 
special  significance  have  been  methods  to  determine  the  IR 
"signatures"  of  military  targets,  especially  aircraft.  The 
Aeronautical  Systems  Division  Infra-Red  Emission  Prediction 
Program  (ASDIR)  is  a  large  computer  model  used  by  ASD  and 
civilian  contractors  in  a  wide  range  of  IR  engineering 
applications.  Using  a  geometric  representation  of  the  sub¬ 
ject  aircraft,  along  with  engine  physical  characteristics 
and  operating  parameters,  the  program  computes  emissions 
from  hot  exhaust  surfaces  and  the  exhaust  plume.  ASDIR 
then  estimates  the  propagation  of  that  energy  through  an 
attenuating  atmosphere  to  an  observer  (usually  an  IR 
detector).  With  specific  detector  characteristics,  a 
prediction  can  be  made  as  to  whether  the  target  could  be 
detected  at  a  given  range,  altitude  and  viewing  angle. 

A  serious  complicating  factor  in  such  predictive 
schemes  is  atmospheric  attenuation  of  the  IR  energy. 
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The  energy  is  absorbed  and  scattered  by  both  the  molecules 
of  the  air  and  by  atmospheric  aerosols  such  as  dust,  smoke, 
and  salt  particles.  The  total  effect  of  the  atmospheric 
constituents  can  be  approximated  by  the  effects  of  the 
separate  attenuators.  The  attenuating  mechanisms  of  molec¬ 
ular  scattering,  molecular  absorption  (especially  by  ^0 
and  CO2) ,  and  aerosol  scattering  and  absorption  can  be 
related  to  meteorological  parameters  such  as  temperature, 
pressure,  relative  humidity,  and  particle  density.  Then, 
in  theory,  one  can  predict  the  cumulative  attenuation . of 
radiation  as  a  function  of  the  path  and  meteorological 
observables . 

Aerosol  attenuation  is  the  least  understood  of  the 
attenuating  mechanisms  listed  above.  The  aerosol  environ¬ 
ment  in  any  location  is  a  result  of  complex  aerosol 
production  and  removal  processes,  relative  humidity,  and 
particle  composition.  Realistic  aerosols  are  difficult  to 
describe  and  model,  and  yet  attenuation  by  aerosols  is 
often  the  limiting  factor  in  IR  detection,  especially  at 
low  altitudes  where  aerosols  are  concentrated.  When  ASDIR 
was  composed  no  provisions  were  made  to  compute  attenuation 
by  aerosols.  While  the  program  gives  good  results  (as 
compared  against  actual  measurements)  at  high  altitudes  and 
high  latitudes,  its  predictions  at  lower  altitudes  are 
unrealistic. 

LOWTRAN,  another  large  computer  program  developed  by 
the  Air  Force  Geophysics  Laboratory,  predicts  atmospheric 
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IR  transmission  over  the  spectral  range  from  0.25  to 
28.5  pm.  LOWTRAN  computes  attenuation  by  all  the  major 
atmospheric  constituents  and  even  by  several  of  the  trace 
gases  in  air.  The  latest  edition  of  LOWTRAN  includes  five 
aerosol  models  and  makes  provision  for  the  user  to  insert 
his  own  aerosol  model  if  desired.  LOWTRAN  is  the  most 
widely  distributed  transmission  code  and  has  become  a 
standard  in  the  field  of  atmospheric  propagation  studies. 

Objectives 

The  objectives  of  this  independent  study  are  to  write 
an  aerosol  subroutine,  based  on  the  physics  of  LOWTRAN,  for 
insertion  into  the  ASDIR  code  and  to  investigate  the  atten¬ 
uation  of  IR  radiation  by  atmospheric  aerosols.  A  minor 
modification  to  the  LOWTRAN  logic  is  made  to  account  for 
results  from  a  recent  attenuation  study  at  Grafenwohr, 
Federal  Republic  of  Germany. 

Approach 

LOWTRAN  parameters  needed  to  predict  aerosol  trans¬ 
mission  are  (1)  assumed  aerosol  particle  densities 
(2)  equivalent  sea  level  path  lengths  for  vertical  paths, 
AV(z)  and  (3)  total  sea  level  extinction  coefficients  for 
the  "average  continental"  and  "maritime"  aerosol  models. 

For  this  study,  curve  fits  were  obtained  for  these  param¬ 
eters  giving  analytical  expressions  which  are  easy  and 
efficient  for  use  in  the  large  ASDIR  computer  code. 

Using  the  curve  fit  results,  a  subroutine  called  HAZE 
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was  written  which  computes  the  aerosol  transmission  as  a 
function  of  wavelength,  meteorological  range,  altitudes, 
and  slant  range  for  either  of  two  selected  aerosol  models. 
The  subroutine  was  validated  against  LOWTRAN  and  installed 
into  ASDIR.  Several  runs  of  the  modified  ASDIR  program 
were  made  using  an  assumed  black  body  "target"  at  a  given 
temperature  and  range.  The  total  transmission  curves  pro¬ 
duced  by  the  modified  ASDIR  were  compared  against  similar 
predictions  from  LOWTRAN  to  validate  the  installation  of 
the  HAZE  subroutine  into  ASDIR. 

Assumptions 

Several  specific  assumptions  that  should  be  mentioned 
here  are: 

(1)  The  aerosol  extinction  coefficients  used  in 
LOWTRAN  result  from  Mie  theory  calculations  which  assume 
the  particles  are  monodispersed,  homogeneous  spherical 
particles.  Most  real  aerosols  consist  of  nonspherical 
particles . 

(2)  LOWTRAN  particle  densities  are  characteristic  of 
only  a  few  real  atmospheres.  Actual  densities  at  any 
location  may  vary  considerably  from  the  LOWTRAN  distribu¬ 
tions  because  of  local  aerosol  production  mechanisms  (eg. 
factories,  desert  sand,  ocean,  etc.) 

(3)  Beer's  Law  is  used  in  this  study  to  compute  the 
transmission  of  energy  in  a  narrow  waveband  (typically 
20  cm-1).  The  law  holds  exactly  only  at  a  single 
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(monochromatic)  wavelength.  Here,  an  average  extinction 
coefficient  over  the  small  waveband  is  used  to  decrease  the 
running  times  of  the  computer  programs .  Higher  accuracy 
can  be  achieved  by  decreasing  the  wavelength  interval  but 
with  a  corresponding  increase  in  cost. 

Development  of  the  Report 

In  Chapter  II  a  brief  description  of  the  ASDIR  emis¬ 
sion  program  is  given  to  further  define  its  purpose, 
methodology,  and  its  need  for  an  aerosol  prediction  capa¬ 
bility.  Next,  a  general  treatment  of  atmospheric  aerosols 
is  given  with  emphasis  on  the  quantities  used  by  the 
LOWTRAN  aerosol  predictive  scheme.  A  brief  section  on  the 
theory  of  IR  extinction  follows  in  which  Beer's  Law  is 
defined  along  with  the  relationship  between  meteorological 
range  and  particulate  density.  Chapter  II  closes  with  a 
brief  discussion  of  LOWTRAN,  again  with  emphasis  on  its 
aerosol  attenuation  prediction  scheme.  The  development  of 
the  subroutine  HAZE  for  inclusion  into  ASDIR  is  described 
in  Chapter  III.  Results  of  the  systematic  curve  fits  of 
extinction  coefficients,  particle  densities,  and  aerosol 
equivalent  path  lengths  are  listed.  Results  of  the  sub¬ 
routine  HAZE  calculations  are  presented  in  plots  and 
compared  against  LOWTRAN  results  for  similar  cases. 

Chapter  III  also  discusses  the  results  obtained  with  the 
modified  ASDIR  program  and  the  value  of  an  aerosol 
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computation  capability.  Comparisons  are  given  between 
modified  ASDIR  and  LOWTRAN  results  for  a  few  cases.  Limi¬ 
tations  of  the  modified  ASDIR  program  are  discussed  along 
with  recommendations  for  future  development  of  the  ASDIR 
code . 
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II.  Background  and  Theory 

ASDIR  -  IR  Emission  Prediction  Program 

The  Aeronautical  Systems  Division's  Infra-Red  Emission 
Prediction  Model  (ASDIR)  was  developed  in  1972-74  and  is 
documented  in  three  volumes  (Stone  and  Tate,  1975a,  b,  c) 
available  through  the  Defense  Documentation  Center  (DDC) . 
ASDIR  is  a  composition  of  programs  that  compute  the  IR 
energy  emitted  by  the  airframe,  engine,  and  plume  of  an 
aircraft  in  flight. 

Infrared  energy  is  constantly  being  emitted  by  all 
objects.  In  a  complex  object  such  as  an  aircraft  in  flight 
important  emitters  are: 

(1)  hot  parts  of  engine  exhaust  systems 

(2)  hot  gases  of  engine  exhaust  (plume) 

(3)  aerodynamical ly  heated  leading  edges  and  surfaces 

(4)  heat  exchangers 

(5)  aircraft  lights 

In  addition,  background  radiation  from  earth,  sky,  clouds, 
etc. ,  and  reflection  of  sunlight  may  be  very  important  at 
certain  viewing  aspects  and  in  certain  wavelength  bands. 
Figure  1  shows  pictorially  the  various  sources  of  IR  radi¬ 
ation  emanating  from  an  aircraft. 

The  energy  radiated  (or  reflected)  from  an  aircraft 
will  consist  of  both  gray-body  continuum  radiation  and 
molecular  species  (line  spectra)  radiation  from  the  plume. 
IR  emissions  of  a  subject  aircraft  are  summed  from  external 
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Aircraft  IR  Radiation  Sources 


surface  and  plume  calculations  over  the  frequency  spectra. 
The  summed  emissions  are  presented  as  aspect  distributions 
in  azimuth  and  elevation  of  IR  energy  emitted  from  the 
aircraft  and  transmitted  through  the  atmosphere  to  an  ob¬ 
server  at  various  ranges.  The  relative  importance  of  the 
various  contributors  to  the  IR  signature  is  shown  in  Fig  2. 
ASDIR's  treatment  of  each  of  the  main  energy  contributors 
and  the  atmospheric  transmission  scheme  will  now  be 
considered. 

Hot  Surface  Emissions.  ASDIR  uses  a  subroutine, 

SIGNIR,  to  predict  IR  emissions  from  axisymmetric  turbojet, 
turbofan,  and  turboshaft  engine  exhaust  system  hot  surfaces. 
SIGNIR,  a  product  of  the  Vought  Aeronautics  Company,  con¬ 
ducts  an  internal  flow  analysis,  includes  special  surface 
cooling  effects,  and  computes  a  thermal  balance  to  provide 
the  system  surface  temperature  distribution  (Stone  and 
Tate,  1975c: 2).  The  routine  then  calculates  surface  radia¬ 
tion  (with  multiple  surface  reflections)  which  emanates 
from  the  nozzle  exit  of  the  engine.  The  radiant  energy  is 
calculated  in  10  to  25  wavenumber  increments  (wavebands) 
for  IR  wavelengths  from  0.9  pm  to  200  ym. 

A  flowchart  for  SIGNIR  is  shown  as  Fig  3.  Inputs  to 
the  routine  include  the  emissive  properties  of  the  metallic 
surfaces  and  the  relative  physical  dimensions  of  engine 
exhaust  system  components.  After  a  one-dimensional  compres¬ 
sible  flow  analysis  is  performed,  boundary  layer  calcula¬ 
tions  are  made  to  obtain  coefficients  of  convection  and 
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Relative  Importance  of  Contributions  to  a  Weapon  Systems  Infrared  Signature 

(Ref  JTCG/AS  Handbook,  1977) 


Figure  3.  SIGNIR  Flowchart  as  Used  in  ASDIR 
(Ref.  Stone  and  Tate,  1975c: 3) 
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friction  along  each  surface  in  the  system.  Next,  any  cool¬ 
ing  effects  are  considered  such  as  transpiration,  film  or 
convection- film  methods.  Internal  geometric  view  factors 
are  then  calculated  which  take  advantage  of  symmetry  and 
account  for  shadowing  when  the  view  between  points  is 
obscured  by  another  surface. 

After  the  geometric  view  factors  are  computed,  and 
with  surface  emissivities  available  as  input  data,  an  over¬ 
all  radiation  exchange  factor  is  calculated.  This  factor, 
used  in  the  radiation  calculations,  accounts  for  the  view 
between  two  surfaces,  the  non-black  nature  of  the  surfaces 
and  the  multiple  reflections  occurring  between  all  surfaces. 
The  procedure  assumes  that  all  the  emitting  surfaces  are 
Lambertian. 

SIGNIR  performs  an  internal  heat  balance  on  the  system 
which  considers  radiation,  convection  and  conduction.  This 
balance  provides  the  temperature  at  each  of  the  internal 
surface  nodes  of  the  exhaust  system  for  the  calculation  of 
the  hot  parts  energy  radiating  from  the  exit  of  the  engine 
exhaust.  The  radiant  energy  from  the  internal  engine  hot 
parts  is  considered  to  represent  the  background  energy  for 
those  "rays"  which  emit  to  the  observer  from  the  engine 
nozzle  exit.  (See  plume  discussion  for  energy  attenuated 
by  the  plume  gases) . 

External  Surface  Emissions.  External  surfaces  of  the 
aircraft  are  treated  as  independent  surfaces  which  occlude 
background  incident  energy  and  present  their  own  grey-body 


12 


Lambertian  energy.  Each  of  an  array  of  up  to  20  surfaces 
is  characterized  by  area,  temperature,  and  emissivity. 

This  array  is  prepared  separately  and  provided  to  the  ASDIR 
program  as  input.  This  treatment  represents  all  external 
aircraft  surface  emitters  such  as  aerodynamically  heated  or 
cooled  regions  (heat  exchangers  and  exhaust  ducts),  lights, 
and  even  sun  glint  spots . 

Plume  Emissions  and  Atmospheric  Attenuation.  These 
two  apparently  independent  quantities  are  discussed  as  if 
they  were  common  because  they  are  treated  together  within 
ASDIR.  A  subroutine  (PLUSIG)  treats  the  gaseous  emission 
and  absorption  of  molecular  species,  accounts  for  line 
broadening  due  to  pressure  and  optical  depth,  and  is 
designed  to  process  these  calculations  over  a  wide  range 
of  temperatures.  The  calculation  procedure  for  PLUSIG  was 
worked  out  between  Dr.  G.  Lindquist  of  The  University  of 
Michigan  and  Dr.  C.  Stone  of  ASD. 

Background  energy  is  treated  as  entering  the  plume 
"ray"  from  a  point  away  from  the  observer  and  contributes 
to  the  energy  after  absorption  in  the  plume  gases.  Should 
the  background  energy  not  be  fully  attenuated  (absorbed) 
in  passing  through  the  plume,  the  plume  appears  as  trans¬ 
parent.  This  transparency  is  most  prominent  for  those 
plume  rays  which  emit  from  the  engine  exhaust  nozzle. 

Plume  rays  are  laid  out  on  a  geometrical  basis  taking 
advantage  of  axisymmetric  symmetry  (See  Fig  4  for  plume  ray 
geometry, .  The  origin  of  the  ray  geometry  is  the  center  of 
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Projection  of 
emitting  element 


Figure  4.  ASDIR  Ray  Geometry 


the  engine  exhaust  nozzle  exit.  Rays  are  divided  into  NT 
segments,  NA  of  which  are  atmospheric.  ASDIR  uses  a  paral¬ 
lel  ray  assumption  (treating  the  aircraft  as  a  point  source 
which  was  a  simplifying  approximation  of  greater  benefit  to 
the  treatment  than  the  loss  of  accuracy  it  invokes.  This 
consideration  was  based  on  included  angles  and  relative 
ranges  of  intended  typical  problems. 

From  this  short  discussion  it  can  be  appreciated  that 
the  calculation  procedure  is  complex  and  lengthy  in  view  of 
the  highly  dependent  nature  of  plume  radiation  on  gaseous 
species  concentration  and  distribution  in  the  plume,  the 
plume  geometry,  and  absorption-emission  physics  of  the 
plume  gases.  The  general  approach  to  the  problem  taken  by 
PLUSIG  is  a  gas  band  model  based  on  the  Curtis-Godson 
approximation.  This  technique  is  described  in  some  detail 
by  Goody  (1964:236ff)  and  involves  assigning  to  a  non- 
homogeneous  path  (such  as  encountered  in  the  plume)  a 
temperature-scaled  amount  of  absorbing  matter  at  a  mean 
pressure.  A  single  spectral  line  absorbing  along  a  non- 
homogeneous  path  can  then  be  simulated  by  a  line  absorbing 
along  a  homogeneous  path  with  both  the  intensity  and  width 
of  the  line  empirically  varied.  Band  parameters  of  line 
intensity,  widths,  and  separations  of  the  lines  are  func¬ 
tions  of  both  wavelength  and  temperature  for  band  models 
which  both  emit  and  absorb.  These  band  parameters  were 
experimentally  determined  in  the  laboratory  from  absorption 
measurements  in  isothermal  samples  of  gas.  Their  use  in  a 
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computer  code  involves  an  array  of  the  band  parameters 
which  are  adjusted  to  approximate  a  homogeneous  segment  of 
gas  in  a  plume.  For  the  (assumed)  homogeneous  segment  of 
gas  the  Planck  emission  law  is  then  used  to  calculate  its 
radiant  energy. 

PLUSIG  divides  the  total  path  through  the  plume  to  the 
observer  into  NT  segments,  NA  of  them  atmospheric  as  shown 
in  Fig  4.  Each  of  the  NT  segments  along  a  ray  is  described 
by  the  gaseous  partial  pressures,  temperatures,  and  length. 
The  transmittance  of  each  ray  segment  is  computed  by  PLUSIG 
using  several  other  subroutines.  In  the  atmospheric  seg¬ 
ments  PLUSIG  includes  only  molecular  I^O,  CC>2 ,  and  ^0  as 
absorbing  gases  excluding  aerosol  effects.  Work  in  this 
independent  study  was  directed  toward  providing  the  capa¬ 
bility  to  include  aerosol  attenuation  in  the  atmospheric 
transmittance  calculations  of  PLUSIG  in  the  ASDIR  program. 

Specifically,  a  subroutine  named  ATMOS  is  addressed  by 
PLUSIG  which  calculates  the  transmittances  (•m's)  for  each 
of  the  NA  atmospheric  segments  along  a  ray  from  the  subject 
aircraft  to  an  observer.  Using  the  methodology  developed 
in  LOWTRAN  and  described  in  Chapter  II  of  this  report,  a 
subroutine  was  written  which  is  called  by  ATMOS  for  modifi¬ 
cation  of  the  atmospheric  transmittances  as  necessary  to 
account  for  aerosol  extinction.  Extinction  represents 
scattered  energy  as  well  as  energy  absorbed  by  the  aerosol 
particles.  The  total  transmittance  is  taken  as  the  simple 
multiple  of  the  value  computed  without  aerosol  extinction 
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and  the  value  computed  by  the  new  subroutine  for  each 
wavenumber- band. 

Atmospheric  Aerosols 

The  atmospheric  air  is  never  free  of  particles  having 
a  wide  variety  of  origins,  sizes,  and  chemical  compositions. 
These  particles  constitute  two  main  categories:  (1)  non- 
hygroscopic  particles  such  as  dust,  volcanic  ash,  and  even 
interplanetary  debris  and  (2)  hygroscopic  particles  that 
readily  absorb  and  combine  with  water  such  as  salt,  hydro¬ 
carbons  exuded  by  foliage,  and  many  products  of  combustion. 
Collectively  these  suspended  particles  are  referred  to  as 
aerosols.  If  we  neglect  clouds,  rain,  and  fog  droplets 
(i.e.  visible  water)  and  consider  only  the  aerosols  sus¬ 
pended  in  "dry"  air,  the  term  "haze"  aerosol  emphasizes  the 
particulate  nature  of  these  efficient  scatterers  of  infra¬ 
red  energy.  Haze  occurs  in  greatest  concentrations  near 
the  earth's  surface  by  virtue  of  relatively  high  molecular 
weight  and  is  usually  the  determining  factor  of  visibility 
(McCartney,  1976:114).  Additionally,  and  of  prime  impor¬ 
tance  in  this  study,  haze  particles  constitute  an  important 
source  of  IR  attenuation  throughout  a  wide  spectrum  of 
wavelengths.  Thus,  mathematical  modeling  of  haze  forma¬ 
tions  and  prediction  of  haze  attenuation  are  critical 
elements  in  estimating  atmospheric  transmission  of  radiant 
energy  through  "real-world"  atmospheres. 

The  following  sections  provide  a  basis  for  the 
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theoretical  treatment  of  haze  attenuation  which  follows. 
Here  are  discussed  sources  of  haze  aerosols  and  particle 
composition,  particle  removal  processes,  particle  size  and 
size  distributions,  the  complex  index  of  refraction  of  haze 
aerosols,  and  the  aerosol  vertical  distributions  developed 
for  use  in  the  LOWTRAN  transmission  mathematical  model. 

Sources  of  Atmospheric  Aerosols.  Particles  in  the 
atmosphere  are  formed  by  a  variety  of  processes  some  of 
which  are  man-made  and  others  strictly  natural  phenomena. 

In  some  cases  it  is  difficult  and  perhaps  rather  arbitrary 
to  try  to  distinguish  between  the  source  types.  Cadle 
(1966:7)  classifies  the  formation  mechanisms  into  five  main 
groups : 

(1)  Condensation  of  vapor  and  sublimation  of  solids 
and  the  formation  of  smoke  by  both  natural  and  man-made 
processes . 

(2)  Chemical  reactions  involving  trace  gases  in  the 
atmosphere . 

(3)  Mechanical  disruptions  and  dispersal  of  matter 
from  the  earth's  surface  (includes  sea  salt  over  the  ocean 
and  various  dusts  over  the  continents). 

(4)  Coagulation  of  fine  particles  to  form  larger 
particles  of  mixed  composition. 

(5)  Influx  of  interplanetary  debris. 

The  diverse  sources  of  haze  aerosols  produce  widely 
differing  sizes  and  concentrations  of  particles  which  are 
highly  dependent  upon  the  local  conditions  (eg.  urban, 


18 


rural,  or  maritime  locales)  and  consequently  extremely 
complex  to  empirically  represent.  An  additional  complica¬ 
tion  is  induced  by  the  great  variation  in  the  attenuation 
effects  from  hygroscopic  to  non-hygroscopic  particles.  The 
models  generally  encountered  in  the  literature  usually 
empirically  describe  the  aerosol  concentrations  as  a  speci¬ 
fied  sum  of  percentages  of  dry  and  hygroscopic  particles. 

Dust-Like  (Non-Hygroscopic)  Particles.  Dust-like 
particles  are  continually  being  placed  in  the  atmosphere  by 
several  processes.  As  the  earth  sweeps  through  space, 
interplanetary  debris  that  range  in  size  from  small  meteor¬ 
ites  up  to  several  tons  may  partially  disintegrate  upon 
entering  the  atmosphere  and  shower  the  earth's  surface  with 
millions  of  small  particles. 

Volcanoes  inject  ash  and  soot  into  the  atmosphere, 
sometimes  with  spectacular  results.  McCartney  (1976:116) 
comments  that  the  eruption  of  Krakatoa  in  1883  raided 
clouds  of  ash  to  90,000  ft  and  completely  darkened  the 
skies  for  100  miles.  Within  a  few  weeks  the  smaller  par¬ 
ticles  from  the  eruption  had  completely  circled  the  globe 
at  high  altitudes  and  produced  unusual  visual  effects  at 
sunrise  and  sunset  for  several  years. 

Soil  dust  blown  into  the  atmosphere  by  the  sculpting 
and  lifting  forces  of  winds  of  even  moderate  speed  may 
often  be  transported  great  distances  by  atmospheric  winds 
aloft.  Forest  fires  are  prolific  manufacturers  of  aerosol 
haze  particles  some  of  which  are  non-hygroscopic. 
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Cadle  (1966:12)  estimated  that  an  ordinary  grass  fire 
covering  only  one  acre  yields  2x10 2 2  fine  particles  ranging 
in  composition  from  inorganic  ash  to  complex  tars  and  resin. 
If  all  of  those  particles  were  distributed  uniformly  over 
the  acre  to  an  altitude  of  10,000  ft  the  resulting  concen¬ 
tration  would  be  about  2x10 9  cm"3.  As  in  the  case  of 
volcanic  ash  and  desert  dust,  the  residue  of  fires  may 
travel  great  distances  in  the  atmosphere. 

Of  ever-increasing  importance  is  the  injection  into 
the  air  of  many  types  of  industrial  products.  Grating, 
blasting,  smelting  and  other  operations  produce  several 
forms  of  dust-like  particles  which  vary  greatly  in  size  and 
number  density.  The  larger  particles  quickly  settle  back 
to  earth  while  the  smaller  (radius<1.0  pm)  move  with  winds 
and  remain  suspended  until  removed  by  one  of  the  processes 
of  atmospheric  scrubbing  described  in  a  later  section. 

Hygroscopic  Particles .  The  most  significant  type  of 
hygroscopic  particle  in  the  atmosphere  is  sea  salt  composed 
mainly  of  sodium  chloride.  The  importance  of  sea  salt 
particles  results  from  their  being  so  numerous;  they  are 
widespread  over  both  oceans  and  continental  areas.  Besides 
being  a  consideration  in  radiant  energy  attenuation,  sea 
salt  particles  also  play  a  role  as  nuclei  in  cloud  forma¬ 
tion.  Byers  (1965:70)  reports  that  a  nucleus  of  NaCl  of 
10" 9  gm  dry  mass  (approximately  5  pm  equivalent  spherical 
radius)  becomes  a  water  drop  of  25  pm  radius  at  a  relative 
humidity  of  99  per  cent.  Thus,  sea  salt  particles  in  the 
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atmosphere  present  extreme  complexities  in  aerosol  modeling 
by  their  growth  capability  and  subsequent  changes  in  atten¬ 
uation  effects. 

It  has  been  established,  as  reported  by  McCartney 
(1976:119),  Cadle  (1966:8)  and  others,  that  sea  salt  parti¬ 
cles  are  mainly  formed  by  the  bursting  of  air  bubbles  as 
they  reach  the  surface  of  the  sea.  These  droplets,  after 
evaporation,  produce  particles  with  an  effective  radius  of 
0.15-1.15  Mm.  The  smaller  salt  particles  are  then  caught 
up  by  the  wind,  carried  aloft  and  can  be  transported  great 
distances.  This  explains  why  several  researchers  estimate 
the  aerosol  composition  in  central  Europe  (500  miles  inland) 
to  be  70-80  per  cent  maritime  (largely  sea  salt)  and  20-30 
per  cent  continental  (dust-like)  particles. 

There  are  many  man-made  sources  of  hygroscopic  aerosol 
particles.  Several  types  of  particles  are  produced  by 
photochemical  reactions  and  nucleations  between  combustion 
gases  and  atmospheric  trace  gases.  For  example,  sulfur 
dioxide,  a  result  of  the  combustion  of  coal  and  petroleum 
derivatives,  is  oxidized  by  solar  photosynthesis  to  form 
sulfuric  trioxide.  These  particles  are  very  hygroscopic, 
form  sulfuric  acid,  and  can  become  extremely  corrosive. 

Other  materials  such  as  nitrogen  dioxide  (from  automobile 
exhausts),  hydrogen  sulfide  and  other  sulfate  particles 
abound  in  the  atmosphere. 

Particle  Removal  Processes.  The  particles  being 
placed  in  the  atmosphere  from  the  sources  previously 


21 


described  are  continually  being  removed  or  scrubbed  by 
several  removal  processes.  These  must  be  briefly  discussed 
so  that  an  understanding  of  the  dynamic  nature  of  atmos¬ 
pheric  aerosol  content  may  be  achieved. 

The  simplest  removal  process  is  that  of  gravity  which 
causes  fall-out  or  sedimentation  of  the  particles.  This 
process  is  quite  important  since  it  controls  somewhat  the 
upper  limit  of  aerosol  particle  sizes  (Junge,  1962:123). 

In  settling  through  still  air  a  particle  will  obtain  a  ter¬ 
minal  velocity  when  the  gravitational  force  is  balanced  by 
viscous  drag.  The  viscous  drag  is  a  function  of  air 
density,  so  that  a  particle  having  an  appreciable  fall 
speed  at  a  height  of  80  km  may  essentially  float  at  20  km. 
Sedimentation  may  sometimes  result  in  concentrations  in 
specific  layers,  especially  when  the  particles  are  out  of 
reach  of  rain  which  causes  "washout". 

Washout  is  the  removal  of  particles  by  precipitation 
and  is  very  effective  in  clearing  the  air.  Washout  implies 
some  kind  of  interaction  (joining)  between  the  precipita¬ 
tion  droplet  and  the  aerosol  particle.  Actual  physical 
contact  between  the  two  may  not  be  necessary  and,  in  fact, 
aerosol  particles  are  known  to  be  carried  around  raindrops 
in  the  associated  airstream  (Byers,  1965:76). 

Coagulation  is  a  third  removal  process  and  occurs  when 
two  particles  collide  and  coalesce,  resulting  in  fewer  but 
larger  particles.  The  general  theory  for  coagulation  is 
founded  on  a  statement  of  the  rate  of  collisions  of 
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particles  as  they  diffuse  through  the  air  by  Brownian  or 
turbulent  motion.  An  interesting  result  of  the  theory  (see 
Junge,  1963  for  a  detailed  treatment)  is  that  the  concen¬ 
tration  of  smaller  (Aitken)  particles  is  decreased  very 
rapidly  by  coagulation  as  compared  to  the  larger  particles. 
This  fact  helps  to  fix  the  practical  lower  limit  of  aerosol 
particle  sizes  in  "well-developed"  atmospheres.  A  more 
detailed  treatment  of  aerosol  sizes  and  distributions  is 
presented  in  a  later  section  of  this  chapter. 

Particle  Composition.  It  is  apparent  from  the  pre¬ 
vious  discussions  that  a  wide  diversity  exists  in  sources 
of  aerosol  particles  and  thus  in  particle  composition.  It 
is  therefore  an  extremely  difficult  task  to  estimate  the 
particle  composition,  size  distribution,  and  number  density 
of  aerosols  in  any  particular  application.  On  the  other 
hand,  aerosol  chemical  composition  is  very  important  in 
determining  the  complex  refractive  index  for  use  in  Mie 
theory  calculations  of  scattering  and  absorption  coeffi¬ 
cients.  Many  investigations  have  sought  to  determine  the 
compositions  of  aerosols  and  the  combined  effects  of  the 
production  and  removal  processes.  These  investigations 
have  been  seriously  hampered  by  the  difficult  task  of  meas¬ 
uring  the  composition  of  the  very  small  (Aitken)  particles. 
While  the  approximate  concentrations  of  Aitken  nuclei  can 
be  determined  by  such  means  as  the  Aitken  counter  (similar 
to  a  Wilson  cloud  chamber)  exact  compositions  of  these 
particles  remain  an  open  question.  However,  observations 
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of  the  aggregate  behavior  of  Aitken  particles  (eg.  growth 
curves)  infer  that  their  chemical  composition  is  similar  to 
large  or  giant  particles  and  mainly  of  continental  origin 
(Junge,  1963:157). 

In  summary,  the  constitution  of  aerosol  particles  in 
the  atmosphere  varies  between  two  extremes:  dry  insoluble 
dust  particles  and  completely  soluble  (hygroscopic)  mate¬ 
rials.  Additionally,  any  specific  particle  is  changing  in 
time  under  the  influence  of  coalescence,  coagulation,  and 
combination  with  atmospheric  gases.  These  facts  lead  to 
the  simplifying  step  of  treating  the  atmosphere  by  using  a 
mixed  particle  concept.  This  concept  assigns  an  average 
composition  and  average  index  of  refraction  to  the  entire 
distribution  of  particles.  The  value  of  this  concept  in 
applications  is  obvious;  it  limits  having  to  determine  an 
explicit  aerosol  distribution  for  each  particle  of  differ¬ 
ing  composition.  The  mixed  particle  concept  is  employed  in 
LOWTRAN  calculations  of  extinction  coefficients  and  thus 
represents  a  primary  assumption  in  the  analysis  adapted  for 
this  thesis. 

Particle  Sizes.  An  understanding  of  the  variety  of 
aerosol  production  and  removal  mechanisms  has  been  estab¬ 
lished,  so  it  is  seen  that,  just  as  the  particle  composi¬ 
tions  vary  greatly,  so  does  the  range  of  particle  sizes. 

For  most  applications  in  IR  engineering,  particles  must  be 
considered  whose  radii  extend  over  nearly  four  orders  of 
magnitude.  Figure  5  adapted  from  Junge  (1963:112)  shows 
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both  the  nomenclature  and  importance  of  particle  sizes  for 
the  various  fields  of  atmospheric  science.  Of  special 
interest  in  this  thesis  are: 

Aitken  particles  0.02<r<0.1  ym 

Large  particles  0.10<r<1.0  ym 

Giant  particles  1.00<r<100  ym 

The  various  cutoffs  in  the  particle  radii  are  not  arbitrary 
but  are  based  on  many  measurements  which  yielded  surpris¬ 
ingly  similar  results. 

The  primary  factor  that  determines  the  lower  limit  of 
particle  radii  is  the  rate  of  coagulation  of  the  very 
smallest  particles  (Byers,  1965:72).  Junge  (1963:130)  cal¬ 

culated  the  effects  of  coagulation  and  "aging"  on  a  sample 
haze  aerosol.  Figure  6  shows  how  the  concentration  of  the 
small  Aitken  particles  decreases  very  rapidly  compared  to 
the  time  scale  of  other  meteorological  processes.  The 
figure  shows  that  in  natural  aerosols  the  effects  of  coagu¬ 
lation  are  such  that  the  maximum  of  the  distribution  falls 
between  0.01  and  0.1  ym.  This  occurs  because  particles 
smaller  than  0.005  ym  can  exist  only  briefly  causing  any 
maximum  below  0.01  ym  to  rapidly  disappear. 

On  the  large  end  of  the  particle  spectrum,  sedimenta¬ 
tion  due  to  gravity  is  the  controlling  factor  in  limiting 
the  upper  radius.  Junge  (1963:123ff)  used  general  diffu¬ 
sion  theory  to  predict  that  particles  of  radius  greater 
than  20  ym  would  quickly  settle  back  to  earth  under  normal 
conditions.  He  also  summarized  results  of  several 
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Figure  6.  Size  Distribution  of  Aerosol  Particles  with 
Coagulation  Affecting  the  Lower  Size  Range 

(Ref  Junge,  1963) 
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investigations  which  showed  that  indeed  the  measured  maxi¬ 
mum  cutoff  was  consistently  between  10-20  ym.  Even  over 
the  oceans,  where  particles  larger  than  20  ym  radius  are 
continually  injected  into  the  air  by  bursting  air  bubbles, 
the  20  ym  cutoff  was  still  observed.  LOWTRAN's  use  of 
0.1  ym  and  100  ym  as  lower  and  upper  particle  size  limits 
appears  well  justified  by  experience. 

Within  the  well-established  limits  of  particle  size 
radii  it  is  possible  to  make  some  general  comments  concern¬ 
ing  the  sizes  of  particles  of  different  compositions. 

Dusts,  smokes,  and  dry  haze  which  produce  the  principal 
optical  scattering  in  the  atmosphere  usually  have  radii 
between  0.1  and  1.0  ym  (Byers,  1965:70).  Scattering  by 
these  small  particles  helps  account  for  the  blue  tint  of  an 
otherwise  clear  sky.  Sizes  up  to  100  ym  or  greater  may 
occur  in  dust  storms  or  in  areas  of  industrial  pollution 
but,  as  noted  before,  these  larger  particles  rapidly  settle 
out  of  the  air.  It  has  also  been  pointed  out  that  hygro¬ 
scopic  particles,  especially  sea  salt,  can  grow  very  rapidly 
with  increases  in  relative  humidity.  Cadle  (1966:9)  points 
out  that  below  about  75  per  cent  relative  humidity  sea 
salt/water  droplets  completely  evaporate  and  leave  particles 
about  1.25  ym  in  radius.  The  significance  of  that  size  is 
that  these  particles  are  small  enough  to  be  carried  aloft 
by  atmospheric  winds  and  penetrate  far  inland.  Thus,  many 
locations  in  central  Europe  may  be  characterized  by  the 
continental/maritime  particle  mixes  as  mentioned  earlier. 
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Measured  Particle  Size  Distributions.  An  important 
characterization  of  any  aerosol  is  its  size  distribution 
function  which  defines  the  manner  in  which  the  particle 
population  is  spread  over  the  range  of  sizes.  While  it  was 
once  thought  that  natural  aerosols  consisted  of  fairly  dis¬ 
crete  groups  of  particles,  recent  improved  instrumentation 
has  shown  that  aerosols  have  essentially  a  continuous  size 
distribution  with  well-established  upper  and  lower  limits 
of  particle  size.  Thus,  particle  measurement  investiga¬ 
tions  usually  produce  either  histograms  or  continuous 
distribution  curves  in  which  the  number  of  particles  within 
a  radius  interval  (Ar  or  dr)  is  represented  by  the  area 
over  that  interval.  Because  of  the  wide  range  of  particle 
sizes  and  concentrations,  logarithmic  scales  are  desired. 
The  techniques  used  to  count  and  measure  particles  provide 
data  on  the  number  of  particles  per  specified  radius  inter¬ 
val  (typically  0.1  or  0.2  ym  in  width).  For  very  narrow 
size  classes  (Ar  small)  a  continuous  distribution  may  often 
be  approximated. 

For  a  continuous  distribution  the  number  of  particles 
per  unit  interval  of  radius  and  per  unit  volume  is 


n(r) 


dN 

dr 


(1) 


where  the  differential  quantity  dN  expresses  the  number  of 
particles  having  a  radius  between  r  and  r+dr,  per  unit 
volume  according  to  the  distribution  function  n(r). 


Formally,  the  total  number  of  particles  per  unit 
volume  is  obtained  by  the  integration: 


N  =  /  n(r)dr  (2) 

o 

In  practice,  cutoff  or  limiting  values  of  radii  are  used  as 
limits  instead  of  zero  and  infinity.  Also,  a  cumulative 
concentration  Nc  may  be  defined  as  the  number  of  particles 
with  radii  less  than  r  by  the  equation 

Nc  =  f  n(r)dr  (3) 


where  rQ  is  the  lower  cutoff  for  particle  radius.  A  par¬ 
ticle  distribution  of  this  type  is  the  basis  of  the  aerosol 
attenuation  calculations  in  LOWTRAN  and  is  adapted  for  use 
in  the  modification  to  the  ASDIR  computer  program. 

A  derivative  of  Eq  2  is  often  used  in  practice  where 
a  log  radius  number  distribution  is  defined  because  of  the 
wide  range  in  particle  sizes  and  concentrations.  The  log 
radius  number  distribution  may  be  written  as 


nL(r) 


dN 

d(log  r) 


( cm"  3 ) 


(4) 


where  N  is  again  the  total  concentration  (particles-cm" 3) 
of  aerosol  particles  with  radii  smaller  than  r.  Attempts 
to  describe  particle  size  distributions  have  led  to  several 


proposed  analytical  expressions  to  represent  the  aerosol 
environment 

A  conmonly-used  distribution,  first  developed  by  Junge 
(1963)  from  experimental  work  performed  on  continental 
aerosols,  may  be  written  as  a  "power-law".  Junge  concluded 
that  the  particle  size  distribution  for  an  "average  conti¬ 
nental  aerosol"  can  be  represented  by  an  exponential  curve 
(see  Fig  7)  of  the  form 


dN 


nL(r)  =  a/T-  -  V 
L  d(log  r) 


=  cr 


(5) 


between  r=0.1  pm  and  r=20  jm  .  In  Eq  5 ,  c  is  an  empir¬ 
ical  constant  whose  value  depends  on  the  specific  concen¬ 
tration  and  b  is  dependent  on  the  distribution.  While  b  is 
variable  it  has  been  found  to  be  approximately  3  for  most 
naturally-occurring  aerosols  (McCartney,  1976:141). 

Equation  5  can  be  put  in  nonlogarithmic  form  by  noting 

that 


d(log  r)  =  0.434  d(ln  r)  and  d(ln  r)  ■  ~  (6) 

Thus,  we  can  write 

n(r)  =  -  0.434  cr"(b+1)  (7) 

dr 

Of  course,  actual  particle  size  distributions  may  vary 
considerably  on  an  individual  basis  from  the  power- law 
approximation.  However,  distributions  of  the  form  given  in 
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dN/d(log  r)  (cm-3) 


Radius  (microns) 


Figure  7.  Average  Size  Distribution  for  Continental 
Aerosols  Near  the  Earth's  Surface 

(Ref  Junge,  1963) 


Eq  7  have  been  found  by  many  independent  investigators  to 
be  a  good  representation  of  aerosols  having  a  wide  variety 
of  origins  and  compositions  (McCartney,  1976:141).  Both 
measured  and  deduced  aerosol  concentrations  have  been 
shown  to  agree  quite  well  with  the  predicted  power-law 
expressions . 

Aerosol  Vertical  Distributions.  For  general  applica¬ 
tions,  such  as  those  for  which  ASDIR  is  used,  one  is 
interested  in  IR  attenuation  along  slant  paths  between  an 
emitter  at  one  altitude  and  a  detector  at  another.  Thus, 
it  is  necessary  to  know  how  aerosol  particle  concentrations 
vary  with  altitude.  This  "vertical  aerosol  distribution" 
is  affected  by  a  number  of  factors,  the  most  important  of 
which  is  normally  the  size  distribution  at  the  surface.  The 
local  aerosol  production  mechanisms  and  the  origin  and  his¬ 
tory  of  the  local  air  mass  are  also  significant  factors. 
Qualitatively,  air  layers  over  the  land  masses  are  the  most 
polluted  (naturally  and  artificially)  in  terms  of  aerosol 
particle  concentrations.  Most  of  this  pollution  is  con¬ 
tained  in  a  "mixing  layer"  that  extends  upward  to  3-5  km 
and  can  be  significantly  affected  by  temperature  and  pres¬ 
sure  profiles  in  the  air  layers. 

The  concept  of  the  mixing  (or  boundary)  layer  is  an 
important  one  in  aerosol  attenuation  studies.  The  lowest 
levels  of  the  atmosphere  are  well  stirred  by  winds  and 
convective  currents,  and  the  mixing  ratio  of  aerosol  parti¬ 
cles  is  fairly  constant  in  these  layers.  The  mixing  ratio 
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may  be  expressed  either  as  the  ratio  of  the  total  mass  of 
particles  to  a  unit  mass  of  air,  or  as  particle  to  molecu¬ 
lar  number  density.  Since  the  density  of  air  decreases 
exponentially  with  altitude,  the  particle  concentration 
should  also  decrease  approximately  exponentially  with  alti¬ 
tude.  Penndorf  (1954)  analyzed  data  from  many  airborne 
measurements  of  the  attenuation  of  solar  flux  and  found 
that  particle  concentration  showed  the  expected  exponential 
decrease  with  altitude  in  the  first  5  km  above  the  surface. 
The  LOWTRAN  assumed  particle  distributions  also  show  this 
exponential  decrease  in  the  lower  5  km.  LOWTRAN 's  particle 
distributions  will  be  adapted  without  change  for  the  modi¬ 
fication  to  ASDIR. 

For  an  exponentially  decreasing  particle  distribution 
the  variance  of  particle  concentration  with  altitude  is 

Nz  =  Nq  exp(-Z/Hp)  (8) 

where 

Nq  =  particle  concentration  (cm-3)  at  ground  level 

H  =  scale  height  for  the  aerosol  (ie.  the  altitude 
P  at  which  the  concentration  is  1/e  times  the 
ground  level) 

N  =  concentration  (cm-3)  at  an  altitude  z 

Penndorf 's  data  showed  Hp  to  be  between  1  and  1.4  km  with 
1 . 2  km  being  a  good  average . 

Using  Eq  8  one  can  define  an  equivalent  sea  level  path 
length  through  the  aerosol  environment  completely  analogous 


to  the  concept  of  equivalent  paths  for  molecular  absorption. 
The  equivalent  path  length  is  simply  the  horizontal  dis¬ 
tance  (at  ground  level)  that  has  the  same  attenuation  as 
the  true  (slant)  path  length.  The  concept  lets  one  evalu¬ 
ate  the  amount  of  particles  and  the  resulting  attenuation 
in  terms  of  a  path  length  at  sea  level  where  the  concentra¬ 
tion  is  either  known  or  can  be  estimated  with  reasonable 
accuracy.  Figure  8  illustrates  the  geometry  of  slant  paths 
in  the  plane  parallel  (flat  earth)  approximation.  Using  R 
to  denote  the  equivalent  path  length: 

R  =  HpSecS (exp(-R^cos6/Hp)-exp(-R2COs0/Hp))  (9-a) 

-  Hp (R/Z2~Z1) (exp (-Z/Hp) -exp (-Z2/Hp) )  (9-b) 

where 

=  altitude  of  point  #1 

Z2  =  altitude  of  point  #2 

R  =  slant  range  between  points  #1  and  #2 

R^  =  slant  range  from  sea  level  to  point  #1 

R2  =  slant  range  from  sea  level  to  point  #2 

6  =  zenith  angle  measured  absolute  from  the  local 

vertical 

Hp  =  scale  height  for  aerosol  as  in  Eq  8 

A  word  of  caution  is  in  order  here  because  the  above 
expression  is  developed  using  a  "flat  earth"  approximation 
with  a  nominal  restriction  on  zenith  angles  of  about  75  deg. 
But  at  large  zenith  angles  transmission  is  poor  through  any 
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haze,  and  large  geometric  ranges  are  seldom  involved.  Thus 
a  cautious  use  of  Eqs  9  may  be  justified  even  for  large 
angles . 

With  the  factors  discussed  in  this  section  taken  into 
account,  many  workers  have  devised  models  of  the  vertical 
profiles  of  hazes  found  in  real  atmospheres.  Using  blends 
of  measured  particle  concentrations ,  attenuation  measure¬ 
ments  and  scattering  theory,  Elterman  (1964)  set  up  such 
a  profile  that  has  been  widely  used  in  subsequent  modeling 
efforts.  His  profile  was  later  revised  (Elterman,  1968) 
to  account  for  long-term  stratospheric  effects  on  aerosol 
distributions  that  were  determined  by  new  measurement 
techniques.  Both  models  are  shown  in  Fig  9.  Note  the 
exponential  decrease  in  concentration  from  200  cm- 3  at  sea 
level  to  about  0.02  cm-3  at  10  km  for  the  1964  model.  A 
modification  of  Elterman' s  1968  model  was  incorporated  into 
LOWTRAN  and  is  further  described  in  a  later  section  describ 
ing  the  LOWTRAN  code . 

Complex  Index  of  Refraction.  Since  particles  in  the 
atmosphere  have  such  a  wide  variety  of  chemical  composi¬ 
tions,  it  is  extremely  difficult  to  determine  the  complex 
index  of  refraction  of  a  real  aerosol.  The  index  of 
refraction  is  a  complex  number  whose  real  part  corresponds 
to  scattering  and  whose  imaginary  part  is  used  to  compute 
absorption  in  Mie  theory  calculations.  It  is  apparent  that 
hygroscopic  particles  capable  of  rapid  growth  by  combining 


Particle  Concentration  (cm-3) 


Figure  9.  Vertical  Profiles  of  Particle  Concentration 
from  Elterman  (1964,  1968) 
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with  water  add  to  the  complexity  of  determining  the  total 
index  of  refraction. 

The  real  part  of  the  index  of  refraction  of  common 
aerosols  has  been  estimated  at  about  1.5  by  several  inves¬ 
tigators,  (Volz,  1972;  Hanel,  1976)  while  the  imaginary 
part,  is  less  well-established.  Fortunately,  absorption  by 
aerosols  is  usually  small  and  the  real  part  of  the  index  of 
refraction  is  usually  adequate  to  predict  attenuation  of 
the  particles. 

Volz  (1972)  computed  tables  of  average  refractive 
indices  based  on  studies  of  typical  distillates  of  rain¬ 
water.  These  values  were  used  for  computing  aerosol  extinc¬ 
tion  coefficients  using  Mie  theory  (Appendix  A) .  The 
resulting  extinction  coefficients  were  used  in  the  first 
LOWTRAN  program.  The  original  aerosol  model  in  LOWTRAN  was 
known  as  the  "average  continental  model"  and  assumed  the 
aerosols  to  be  composed  of  70  per  cent  water-soluble  sub¬ 
stances  and  30  per  cent  dust-like  substances.  This  model 
was  considered  typical  of  Central  European  environments. 
Later  editions  have  included  other  model  aerosols  (eg. 
urban,  rural,  maritime)  which  differ  primarily  in  the 
assumed  mix  of  hygroscopic  and  non-hygroscopic  particulates. 
The  values  tabulated  by  Volz  and  the  extinction  coeffi¬ 
cients  computed  from  them  as  used  in  LOWTRAN  are  adopted 
without  change  for  use  in  this  thesis. 
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Extinction  of  IR  Energy  in  the  Atmosphere 


Atmospheric  attenuation  or  extinction  are  terms  used 
to  express  the  processes  by  which  electromagnetic  energy  is 
scattered  and  absorbed  by  the  atmospheric  constituents  as 
that  energy  travels  from  an  emitter  to  a  sensor.  The 
primary  scatterers  and  absorbers  in  the  atmosphere  are  the 
molecules  of  the  air  itself,  water  droplets  in  the  form  cf 
fog  and  precipitation,  and  aerosols  such  as  particulate 
haze.  The  fundamental  premise  of  attenuation  of  energy  is 
the  Beer-Lambert  law  which  states  that  the  attenuation  is 
dependent  on  the  amount  of  matter  along  the  path  of  the 
radiation.  The  law  assumes  that  the  physical  state  (temper¬ 
ature,  pressure,  composition)  of  the  medium  is  uniform 
along  the  patn.  The  monochromatic  spectral  transmission  of 
optical  radiation  is  then 


where 


t  =  I/IQ  =  exp (-8tL) 


=  transmittance  in  per  cent  of  original  radia¬ 
tion  propagated 

■  incident  intensity  of  radiation 
=  emerging  intensity  of  radiation 
=  total  extinction  coefficient 
=  length  of  path  traversed  by  the  radiation 


The  total  extinction  coefficient  0^  reflects  both 
processes  of  scattering  and  absorption  and  may  be  written 
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(11) 


BSc  + 


Ab 


where 

gSc  =  total  scattering  coefficient 
=  total  absorption  coefficient 


The  two  coefficients  above  (6gc  and  3^)  may  be  further 
broken  down  into  aerosol  and  molecular  scattering  and 
absorption  coefficients  as 


Sc  ^Sc/aero  +  ^Sc/mol 

(12-a) 

Ab  *  ^Ab/aero  +  ^Ab/mol 

(12-b) 

and  the  total  aerosol  extinction  coefficient  may  be  written 
as 


®aero  ^Sc/aero  +  ^Ab/aero 


(13) 


The  values  of  3aero  as  a  function  of  wavelength  of 
incident  radiation  may  be  calculated  by  Mie  theory  as 
briefly  outlined  in  Appendix  A  where  3aero  is  dependent 
upon  the  aerosol  particle  radius,  the  complex  index  of 
refraction  of  the  particle,  and  the  particle  size  distri¬ 
bution.  As  developed  in  Appendix  A,  the  single  most 
important  factor  in  determining  the  aerosol  extinction 
coefficient  is  the  dimensionless  size  parameter  a  given  by 


a  - 


2r 

X 


(14) 
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where 

r  =  particle  radius 

X  =  wavelength  of  incident  radiation 

Hie  calculations  show  that  absorption  by  most  dust-like 
particulates  is  very  small  and  can  be  reasonably  neglected 
in  most  applications.  Scattering  by  all  aerosol  particles, 
however,  is  a  very  important  part  of  atmospheric  attenua¬ 
tion  and  must  be  accounted  for  in  some  detail. 

Considering  the  other  factors  in  Eqs  12,  scattering  by 
the  molecules  of  the  air  itself  is  significant  at  optical 
wavelengths  and  helps  account  for  the  blue  tint  of  the  sky. 
At  wavelengths  greater  than  0.6  ym  (which  includes  the  IR 
portion  of  the  spectrum)  this  scattering  by  the  extremely 
small  molecules  becomes  negligible  compared  to  aerosol 
scattering  (Rensch,  1968:1).  Thus,  molecular  (also  termed 
Rayleigh)  scattering  is  safely  omitted  in  the  ASDIR  program 
for  the  long  wavelength  applications. 

Molecular  absorption,  denoted  by  8^, /moi  in  Eq  12-b 
is  an  extremely  complicated  function  that  represents  a 
major  hurdle  in  computer  prediction  codes  such  as  ASDIR  and 
LOWTRAN.  Both  codes  use  a  "waveband"  model  to  calculate 
molecular  absorption  in  very  narrow  wavelength  bands  for 
the  major  molecular  attenuators.  These  attenuators  in 
order  of  importance  are:  1^0,  CC^,  0^,  ^0,  CO,  O2,  CH^, 
and  N2  (McClatchey,  et  al,  1972:11).  1^0,  CO2,  and  O3 

(in  some  narrow  wavelength  bands)  are  by  far  the  most 
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important  absorbers.  Both  LOWTRAN  and  ASDIR  have  demon¬ 
strated  ability  to  calculate  the  absorption  by  molecular 
gases  to  a  reasonable  accuracy  in  most  areas  of  the 
spectrum.  Calculation  of  aerosol  extinction,  especially  in 
the  lowest  layers  of  the  atmosphere,  is  still  a  large  uncer¬ 
tainty  in  LOWTRAN  and  has  been,  as  noted,  nonexistent  in 
ASDIR. 

Aerosol  Extinction.  The  discussion  thus  far  of  the 
physical  production  mechanisms  of  aerosols  qualitatively 
suggested  the  importance  of  aerosols  as  attenuators  of 
electromagnetic  energy.  This  section  will  treat  the  total 
extinction  due  to  aerosols  as  predicted  by  Mie  theory.  The 
essence  of  this  section  is  the  use  of  the  extinction  coef¬ 
ficients  predicted  by  Mie  theory  and  the  connection  between 
certain  meteorological  parameters  and  aerosol  extinction. 

Mie  theory  predicts  an  extinction  coefficient  for 
spherical  polydispersed  particles  as  (Rensch,  1964:4) 

r2 

6ext  =  /  Qext(m,r)n(r)rzdr  (15) 

rl 

where 

n(r)  =  the  size  distribution  of  the  spheres 

Qoirt- (m> r)  =  extinction  efficiency  at  radius  r 
exC  (See  Eq  45) 

r  =  particle  radius 

m  =  complex  index  of  refraction  of  the  particle 


rl'  r2 


lower  and  upper  limits  of  the  aerosol 
distribution 


Assuming  a  size  distribution  ns  suggested  by  Junge  and 
given  by  Eq  7,  substitution  into  Eq  15  gives 


l'aero  "  °-43Aci'  /  Qcxt  (>’>.  r)  r1_bdr 
rl 


(16) 


Changing  the  independent  variable  to  a*27ir/A  converts 
Eq  16  to  the  following: 


aero 


0 . 4  34  cn 


r 2  7i ' 

b-2 

A 

^xt 

b-1 


dot 


(17) 


where 


0 . 4 34 c  n 


b-2 


'ext 

b-1 


dot 


al 


(18) 


09) 


The  role  of  the  constant  c  in  Kq  17  needs  closer  exam¬ 
ination.  From  Eq  7  it  is  seen  that  the  dimensions  of  c  are 
fa—  ^ 

(length)  *  and  the  value  of  c  depends  on  the  particle 
concentration.  Rearranging  Eq  17  shows  that 


c  -  2.304 


A  |b~2  |Naero 
u  K 
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(20) 


Next,  a  well-established  relation  between  meteorolo¬ 
gical  range  and  the  extinction  due  to  aerosols  and  air 
molecules  at  X=0.55  ym  is  the  Koschmeider  equation 
(McCartney,  1976:42) 


R 


m 


3.91 


Sc 


(21) 


where  absorption  in  the  visual  wavelengths  is  neglected  and 

T,  the  observer's  visual  contrast  threshold  is  taken  as 

0.02.  The  above  equation  largely  eliminates  the  subjective 

factors  present  in  many  operative  definitions  of  visibility 

Rm  is  obtained  by  specifying  that  distance  at  which  a  large 

black  object  can  be  seen  against  a  day  sky  background  with 

T=0.02  .  The  value  for  T  has  been  validated  by  several 

independent  researchers  as  being  appropriately  average. 

B0_  is  approximately  equal  to  the  total  spectral  extinction 
s  c 

coefficient  at  0.55  ym  which  is  due  primarily  to  scattering 
by  both  aerosols  and  the  molecules  of  the  air.  The  value 
of  A=0.55  ym  is  used  as  the  wavelength  of  greatest  sensi¬ 
tivity  of  the  human  eye.  Equation  21  may  then  be  written 
as 


aero 


(0.55  ym)  = 


3.91 


R  (0.55  ym) 
tn 


Pmol(0-55  yn°  (22) 


where  the  terms  in  parentheses  indicate  that  meteorological 
range  refers  to  a  specific  wavelength  of  0.55  ym. 


Mie  theory  establishes  that  the  aerosol  extinction 
coefficient  Baero  is  directly  proportional  to  the  aerosol 
number  density  N(z)  so  we  can  write 

N(z)  =  +  b(z)  (23) 

“m 

where  a(z)  and  b(z)  are  constants  for  any  given  altitude  z. 
The  above  equation  can  be  used,  as  it  is  in  LOWTRAN,  to 
interpolate/extrapolate  to  visual  ranges  other  than  those 
for  which  specific  values  of  N(z)  are  given  as  in  Table  1 
for  visual  ranges  of  5  km  and  23  km.  A  later  section  in 
this  chapter  discusses  further  the  use  of  this  procedure  in 
LOWTRAN  which  is  adapted  for  use  in  ASDIR. 

LOWTRAN- AFGL  Transmission  Model 

The  most  widely-accepted  computer  model  in  use  for 
predicting  atmospheric  transmittances  is  LOWTRAN,  under 
constant  development  and  improvement  by  the  Air  Force 
Geophysics  Laboratory  (AFGL) .  LOWTRAN  is  an  evolutionary 
product  (LOWTRAN  4  is  the  latest  edition)  that  is  updated 
as  new  data  concerning  molecular  and  aerosol  effects  are 
generated.  The  program  is  "user  oriented"  in  that  it  is 
well  documented,  easy  to  operate,  and  thus  well  suited  for 
field  use.  No  attempt  here  is  made  to  fully  describe  the 
LOWTRAN  methodology,  but  the  interested  reader  is  referred 
to  McClatchey,  et  al  (1972);  Selby  and  McClatchey  (1972, 
1975);  and  Selby,  et  al  (1976,  1978)  for  a  complete 
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or 


background  description  and  documentation  of  the  program. 
Especially  useful  is  Optical  Properties  of  the  Atmosphere 
by  McClatchey,  et  al  (1972)  in  which  the  theoretical  basis 
for  the  LOWTRAN  program  is  fully  outlined  with  exception  of 
the  current  aerosol  treatment.  Here  only  a  brief  descrip¬ 
tion  of  the  program  will  be  given  with  particular  emphasis 
on  the  aerosol  attenuation  method. 

LOWTRAN- A  General  Description.  LOWTRAN  computes  the 
transmittance  of  the  atmosphere  over  a  broad  wavelength 
interval  with  wave  number  spectral  resolution  of  20  cm-1 
(waves  per  cm)  .  It  falls  into  the  general  category  of  wave 
band  model  schemes  and  is  based  on  laboratory  transmittance 
measurements  complemented  by  emission  line  theory  molecular 
constants  compiled  and  presented  by  McClatchey,  et  al 
(1972).  The  program  is  capable  of  predicting  transmittance 
along  arbitrary  slant  paths  in  the  wavelength  region  0.25 
to  28.5  um  (ultraviolet  to  middle  IR  portions  of  the 
spectrum).  The  program  allows  a  choice  of  six  reference 
atmospheres  or  direct  inputs  of  measured  (or  assumed) 
meteorological  parameters.  Thus,  LOWTRAN  is  an  extremely 

I  flexible  tool  which  represents  a  major  accomplishment  as  a 

! 

general  purpose  atmospheric  transmission  code. 

LOWTRAN  calculates  transmission  over  a  selected  path 
by  calculating  the  individual  transmittances  of  separate 

i 

i  contributors  as  though  each  were  acting  independently.  In 

order,  the  elements  included  in  the  LOWTRAN  calculations 
are : 
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(1)  water  vapor  line  absorption  from  350  cm-1  to 
14,500  cm-1 

(2)  "uniformly  mixed  gases"  (CO2,  ^0,  CH^,  CO,  N2  and 
©2)  in  the  intervals  500  to  8060  cm"1  and  12970  to  13190  cm"1 

(3)  nitrogen  continuum  between  2080  and  2740  cm-1 

(4)  water  vapor  continuum  between  670  and  1400  cm-1 
and  between  2380  and  2850  cm-1 

(5)  ozone  in  the  IR  between  575  and  3270  cm-1,  in  the 
visible  between  13000  and  23300  cm-1,  and  in  the  UV  for 
wavelengths  shorter  than  27500  cm-1 

(6)  molecular  scattering  for  wavelengths  shorter  than 
2740  cm" 1 

(7)  nitric  acid  (HNO^)  in  the  bands  850-920  cm-1, 
1275-1350  cm"1,  and  1675-1735  cm*1  (LOWTRAN  4  only) 

(8)  aerosol  scattering  and  absorption  at  all  wave¬ 
lengths 

LOWTRAN  steps  down  the  wavelength  scale  in  5  cm-1 
increments  and  computes  individual  transmittances  t-^(v), 
■^(v),  •••xn(v)  where  v  is  the  frequency  in  wavenumbers 
(cm-1).  Note  that  the  spectral  resolution  of  the  wave 
band  model  in  LOWTRAN  is  20  cm"1  so  that  each  t^(v)  is  an 
average  transmittance  over  20  wavenumbers.  The  predictions 
of  LOWTRAN  thus  approximate  the  results  seen  by  a  low  reso¬ 
lution  spectrometer  (Walsh,  1978:3-24). 

After  the  individual  transmittances  are  computed,  the 
overall  average  transmission  at  the  wavenumber  considered 
(i.e.  the  middle  of  the  20  cm"1  interval)  is  taken  as 
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(24) 


t(v)  *  T1(v)t2(v) ♦ * *Tn(v)  =  nTn(V) 

Walsh  (1978:3-24)  discusses  the  validity  of  the  assumption 
implied  in  Eq  24,  ie.  that  average  transmission  at  v  can 
be  taken  as  the  simple  product  of  the  individual  contrib¬ 
utors  at  the  same  wavelength.  He  concludes  that  the 
relation  might  not  be  valid  where  individual  average  trans¬ 
mit  tances  were  low  and  sharp  structure  was  present  within 
the  scale  of  20  cm-1  for  at  least  two  of  the  contributors. 
Walsh  also  cites  experimental  results  which  show  that  the 
expression  holds  reasonably  well  for  atmospheric  gases  in 
the  overall  transmission  range  of  about  30-70  per  cent. 

LOWTRAN  Calculation  of  Individual  Transmittances . 
LOWTRAN  uses  some  graphical  prediction  schemes  originally 
suggested  by  Altshuler  (1961)  which  assign  a  functional 
form  to  the  average  transmittance  over  a  20  cm-1  interval. 
The  functional  forms  and  their  derivations  are  described  in 
McClatchey  et  al  (1972)  and  are  dependent  on  the  amount  of 
absorber  (eg.  water  content),  pressure,  temperature  and 
wavelength.  The  functional  forms  were  put  into  a  kind  of 
nomogram,  an  example  of  which  is  presented  in  Fig  10.  Such 
charts,  for  each  of  the  absorber  species,  were  then  digi¬ 
tized  to  form  the  basis  of  the  LOWTRAN  computer  program. 

The  aerosol  calculations  are  of  major  interest  in  this 
study  and  will  be  considered  now  in  detail. 

LOWTRAN  Aerosol  Treatment.  The  LOWTRAN  program  uses 
two  vertical  aerosol  profiles  that  describe  a  "clear"  and 
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LOWTRAN  Prediction  Chart  for  Water  Vapor  Transmittance  (4-26  ui 

(Ref  McClatchey ,  et  al,  197 


a  "hazy"  atmosphere  corresponding  to  visibilities  (meteor¬ 
ological  ranges)  of  23  km  and  5  km,  respectively.  The  two 
profiles  are  based  on  the  following  assumptions  (Selby  and 
McClatchey ,  1975:12ff): 

(1)  A  particle  size  distribution  similar  to 
Deirmendj ian ' s  Haze  Model  C  but  with  the  large  particle 
cutoff  extended  to  100  ym. 

(2)  The  particle  size  distribution  is  assumed  to 
remain  constant  with  altitude.  Particle  densities  versus 
altitude  are  listed  in  Table  I. 

(3)  The  variation  of  refractive  index  with  wavelength 
was  obtained  from  measurements  by  Volz  (1972)  who  assumed 
specific  mixes  of  continental  (dust-like)  and  maritime 
aerosols . 

(4)  Aerosol  extinction  coefficients  were  computed 
using  Mie  theory  with  the  assumed  particle  size  distribu¬ 
tion  and  Volz's  refractive  index  values.  Extinction 
coefficients  were  computed  to  correspond  to  several  assumed 
mixes  of  continental  and  maritime  aerosols.  Urban,  rural, 
maritime,  tropospheric,  and  average  continental  aerosol 
models  are  available  in  LOWTRAN  4,  the  latest  edition  of 
the  code.  The  calculated  extinction  coefficients  are 
plotted  in  Figs  11  and  12  for  the  average  continental  and 
maritime  aerosol  models,  respectively. 

LOWTRAN  Aerosol  Interpolation/Extrapolation  Scheme. 

The  relation  between  the  total  extinction  coefficient 
and  the  visibility  (or  meteorological  range)  was  discussed 
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Table  I 


Aerosol  Models.  Vertical  Distributions 
for  a  "Clear"  and  "Hazy"  Atmosphere 

(Ref  McClatchey,  1972) 


PARTICLE  DENSITY  , 

N  (PARTICLES  PER  Ctrl'1) 

■  Altitude 

23  km  Visibility 

5  km  Visibility 

flan) 

Clear 

Hazy 

0 

2.  828E+03 

1.  378E+04 

1 

1.  244E+03 

5.  030E+03 

2 

5.  37  1E+02 

1.844E+03 

3 

2.  25GE+02 

6.  73  lE-*02 

4 

1.  192E+02 

2.  453E  +  02 

5 

8.  987E+01 

8.  987E+01 

6 

6.  337E+01 

6.  337E+01 

7 

5.  890E+01 

5.  r90EH01 

8 

6.069E+01 

6.069E+01 

9 

S.  818E+01 

5.  8  lbE-t  0 1 

10 

5.  675E+01 

5.  675E+01 

11 

5.  317E+01 

5.  317E+01 

12 

5.  585E+01 

5.  585E+01 

13 

5.  15CE+01 

0.  156E+01 

14 

5.046E+01 

5.  048E+01 

15 

4.744E+01 

4. 744E+01 

16 

4. 51 1E+01 

4. 5I1E+01 

17 

4.458E+0I 

4. 458E+01 

18 

4.  314E+0I 

4. 314E+01 

19 

3.634E+01 

3.  634E+01 

20 

2.667E+01 

2.667E+01 

21 

1.  933E+01 

1.  933E+01 

22 

1.455E+01 

1.  455E+01 

23 

1.  I13E+0I 

1.  IJ3E+01 

24 

8.  82GE+00 

8.  826E+00 

25  . 

7.429E+00 

7.  429E+00 

30 

2. 238E+00 

2. 238E+00 

35 

5. 890E-01 

5.  890E-01 

40 

1.  550E-01 

1.  550E-01 

45 

4.  082E-02 

4.  082E-02 

50 

1.  078E-02 

f.  078E-02 

70 

5.  550E-05 

5. 550E-05 

100 

1.  969E-08 

1.  969E-08 
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Figure  11.  LOWTRAN  Average  Continental  Aerosol 

Extinction  Coefficients  (Sea  Level  Values) 
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Figure  12.  LOWTRAN  Maritime  Aerosol  Extinction 
Coefficients  (Sea  Level  Values) 
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earlier.  Middleton  (1952)  reported  the  Koschmeider 
relationship  which  holds  well  at  a  wavelength  of  0.55  pm, 
the  wavelength  at  which  the  human  eye  has  its  greatest 
sensitivity.  Equation  27  defines  the  total  extinction 
coefficient  8t  in  terms  of  visibility  as 


6 


t 


3.91 

VIS 


(27) 


where 


VIS  =  the  horizontal  visibility  or  meteorological 
range  (km) 

T  =  the  threshold  value  of  brightness  contrast  for 
the  eye  (T=0.02  for  daylight  and  large  black 
objects  against  the  sky) 


Eq  27  provides  the  link  between  the  concepts  of  "clear" 

(23  km  visibility)  and  "hazy"  (5  km  visibility)  aerosol 
models  as  used  in  LOWTRAN.  The  following  development  will 
make  this  link  more  apparent. 

From  Eq  27  the  total  extinction  coefficient  may  be 
written  as  the  sum  of  an  aerosol  extinction  coefficient  s 

a 

and  a  molecular  extinction  coefficient  B  .  Thus, 

m 


B 


t 


8+6 

a  m 


3.91 

VIS 


(28) 


and 


e  =  3-91 

a  VIS 


m 


(29) 


From  Mie  theory  calculations  the  aerosol  extinction  coeffi¬ 
cient  is  seen  to  be  directly  proportional  to  the  aerosol 
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number  density  N(z).  We  can  then  write 


N(z)  =  +  b(z)  (30) 

VIS 

where  a(z)  and  b(z)  are  constants  at  any  altitude  z.  In 
Eq  30,  b(z)  is  proportional  to  the  molecular  scattering 
coefficient  at  0.55  ym  and  molecular  absorption  at  0.55  ym 
has  been  assumed  negligible.  The  coefficients  a(z)  and 
b(z)  are  determined  by  linear  interpolation  between  the 
values  of  N(z)  for  5  km  and  23  km  visibilities. 

The  limitation  of  the  LOWTRAN  method  discussed  above 
is  the  assumed  prototypical  aerosol  distributions,  indices 
of  refraction,  and  aerosol  extinction  coefficients  derived 
using  those  assumed  parameters.  No  realistic  alternative 
to  this  semi-empirical  approach  has  yet  been  conclusively 
demonstrated.  However,  recent  measurements  at  Grafenwohr, 
FRG,  (Roberts,  1978)  show  that  for  a  particular  aerosol  and 
spectral  region  the  extinction  is  most  critically  dependent 
upon  the  volume  of  particulate  in  the  atmospheric  path  and 
not  so  much  upon  the  detailed  description  of  the  distribu¬ 
tion  function.  Roberts  recommended  that  LOWTRAN' s  maritime 
aerosol  model  be  used  in  all  low  visibility  situations  - 
a  suggestion  incorporated  in  subroutine  HAZE  written  for 
ASDIR. 
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III.  Results  and  Discussion 


The  approach  taken  to  incorporate  the  LOWTRAN  aerosol 
treatment  into  ASDIR  was  one  which  involved  several  compro¬ 
mises.  The  ASDIR  program  is  a  large  and  long-running 
program  so  that  one  of  the  first  requirements  was  to  keep 
the  added  aerosol  subroutines  as  short  as  possible  without 
substantially  sacrificing  accuracy.  Thus  it  was  decided  to 
curve  fit  several  LOWTRAN  aerosol  parameters  to  obtain 
simple  analytical  expressions  which  reduced  both  the  core- 
storage  requirements  and  running  time  of  the  aerosol 
calculations.  Curve  fits  were  obtained  for  (1)  particle 
densities  as  a  function  of  altitude,  (2)  the  extinction 
coefficients  for  the  average  continental  and  maritime  aero¬ 
sol  models,  and  (3)  for  actual  LOWTRAN  calculated  values  of 
an  important  parameter  -  the  equivalent  sea  level  path  for 
vertical  atmospheric  paths . 

Particle  Densities 

The  assumed  aerosol  vertical  distributions  for  LOWTRAN 
given  in  Table  I  are  used  in  the  calculations  of  equivalent 
path  through  the  aerosol  environment.  The  table  values  for 
the  23  km  visibility  ("clear")  case  and  the  5  km  ("hazy") 
case  are  plotted  in  Fig  13.  Note  that  at  altitudes  above 
the  5  km  mixing  layer  the  assumed  particle  densities  are 
equal  for  all  visibility  ranges.  The  values  were  curve  fit 
using  straight  line  segments  (on  semi-log  paper)  to  obtain 


Table  II 


Derived  Expressions  for  Particle  Density 


z  (km) 

N23(z) 

-  "clear" 

N5(z) 

- 

"hazy" 

0<z53 

2828 

X 

10~. 36z 

13870 

X 

10-.437Z 

3<z-4 

1529 

X 

10-.277Z 

If 

4<zf5 

422 

X 

10” • 137z 

II 

5<z-6 

422 

X 

10”137z 

Same 

as 

N23(z) 

6<z<18 

77 

X 

10" ‘ °139z 

Same 

as 

N23(z) 

18<z<30 

3651 

X 

-j^q-  .  107z 

Same 

as 

N23(z) 

30<z 

Assumed 

1  zero 

Same 

as 

N23(z) 

the  functions  listed  in  Table  II.  The  expressions  for 
N23(z)  and  N^(z)  are  used  to  interpolate  for  the  total 
particle  density  N(z)  using  Eq  30  in  which  the  values  a(z) 
and  b(z)  depend  on  N23(z)  and  N^(z)  at  a  given  altitude  as 
follows : 

a(z)  =  (N5(z)  -  N23(z))/(1/5  -  1/23)  (31-a) 

b(z)  =  (N23(z)/5  -  N5(z)/23)/(1/5  -  1/23)  (31-b) 

These  expressions  were  used  in  HAZE  only  for  computing  the 
equivalent  sea  level  path  length  for  horizontal  transmis¬ 
sion  paths  given  by 


R  -  R(N(z)/No) 


(32) 
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where  R  is  actual  slant  range  as  before.  For  the  general 
case  of  slant  paths,  a  different  approach  was  followed  and 
is  described  in  a  later  section. 

Extinction  Coefficients 

Piecemeal  curve  fits  were  obtained  for  the  average  conti¬ 
nental  and  maritime  extinction  coefficients  plotted  in 
Figs  11  and  12,  respectively.  The  curve  fits  were  partially 
obtained  using  the  IMSL  packages  available  on  the  CDC  6600 
computer  at  Wright- Patterson  AFB,  OH.  The  resulting  analyt¬ 
ical  expressions  are  given  as  Tables  III  and  IV  for  the  two 
selected  aerosol  models.  Maximum  errors  between  the  given 
LOWTRAN  values  and  corresponding  values  computed  using  the 
curve  fit  expressions  were  6  per  cent  for  the  average 
continental  model  and  about  5  per  cent  for  the  maritime 
model.  The  small  errors  that  are  seen  in  the  computer¬ 
generated  transmission  plots  for  LOWTRAN  and  subroutine 
HAZE  are  largely  due  to  these  slight  discrepancies  in  calcu¬ 
lated  extinction  coefficients.  However,  the  benefits  of 
the  curve  fits  were  a  reduction  in  core  storage  requirements 
and  logic  simplification  -  both  desirable  features  for 
inclusion  into  ASDIR. 

Equivalent  Sea  Level  Path  Lengths 

For  transmission  along  a  slant  path  LOWTRAN  first  com¬ 
putes  an  equivalent  sea  level  vertical  path  length  (AV)  for 
the  altitudes  Z^  and  Z2  at  each  end  of  the  actual  path. 

This  vertical  equivalent  path  is  a  measure  of  the  aerosol 
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Table  III 

Curve  Fit  Expressions  for  LOWTRAN  Average  Continental 
Aerosol  Extinction  Coefficients  (Sea  Level  Values) 

(Note:  Values  are  rounded  to  three  significant  figures.) 

Wavelength  -  X  (pm)  Extinction  Coefficient  (km-1) 

0.865X52.0 
2.0  <X<7. 2 

7.2  5X<7 . 9 
7.9  5X58.2 

8.2  <X59.0 
9.0  <X59.2 

9.2  <X5l4. 0 
14.0  <X 


-0.064X  +  0.163 

0.001X2  -  0.015X  +  0.058 

-0.008X  +  0.069 

0.001X  +  0.001 

-0.001X2  +  0.171X  -  0.802 

0.024 

0/611- (1.5  x  10" 4) X+139  exp(-X) 
>0.010 


Table  IV 

Curve  Fit  Expressions  for  LOWTRAN  Maritime  Aerosol 
Extinction  Coefficients  (Sea  Level  Values) 

(Note:  Values  are  rounded  to  three  significant  figures.) 


Wavelength  -  X  (pm) 

Extinction  Coefficient  (knT1) 

0.865X52.5 

-0.027X  +  0.167 

2.5  <X53.5 

0.100 

3.5  <x58.2 

0.002X2  -  0.038X  +  0.204 

8.2  <X<9 . 0 

-0.010X2  -  0.177X  -  0.734 

9.0  5x<11.0 

0.003X2  -  0.061X  +  0.380 

11.05x514.0 

0.004X  -  0.020 

14. 0<X 

0.039 

amounts  in  the  atmosphere  above  the  altitude  z  and  is  given 
by  the  equation 


AV 


^  d(z) 

J  ° 

z 


(33) 


The  quantity  AV  is  a  function  of  both  visibility  and  alti¬ 
tude  through  the  dependence  of  N(z)  on  those  two  variables. 
To  accurately  reproduce  the  LOWTRAN  values  for  AV  in  the 
HAZE  subroutine  LOWTRAN  runs  were  made  to  obtain  actual 
values  for  the  somewhat  lengthy  integrations  required  by 
Eq  33.  The  AV  values  were  plotted  (see  Appendix  A)  and 
curve  fits  were  obtained  as  functions  of  VIS  (visibility) 
and  z  (altitude  in  km) .  The  results  of  the  curve  fits  are 
also  given  in  Appendix  A. 

According  to  the  LOWTRAN  documentation  (McClatchey 
et  al,  1972:35)  the  overall  equivalent  path  length  for  a 
slant  path  is  given  by  Eq  34. 

R  =  (AV(Zl)  -  AV(z2))R/(z2-z1)  (34) 

However,  LOWTRAN  calculations  for  R  did  not  match  exactly 
with  the  expected  results  using  Eq  34.  That  is,  using 
LOWTRAN' s  own  calculated  values  of  AV(z^)  and  AV(z2)  for  a 
given  range  R  and  altitudes  z^  and  z2>  LOWTRAN' s  value  of 
R  differed  as  much  as  9  per  cent  from  that  determined  by 
Eq  34.  The  error  was  found  to  be  dependent  upon  VIS. 
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A  formula  was  devised  to  calculate  this  variance  and,  in 
effect,  a  correction  factor  (CFAC)  was  programmed  into 
the  HAZE  subroutine  to  accurately  reproduce  the  LOWTRAN 
values  for  R.  The  listing  of  subroutine  HAZE  is  given  in 
Appendix  C. 

Aerosol  Transmission  Calculations 

To  determine  the  accuracy  of  the  HAZE  subroutine  in 
duplicating  LOWTRAN' s  aerosol  transmission  calculations 
16  selected  cases  were  run  using  two  aerosol  models  and 
varying  visibilities  and  path  lengths.  Visibilities  were 
chosen  2,  5,  10,  and  23  km  representing  a  variance  from 
very  "hazy"  to  "clear"  atmospheres.  Ranges  of  4  and  10  km 
were  selected  as  representative  of  short  to  medium  range  IR 
missile  detector  operating  ranges.  Figures  14-29  give  the 
plotted  results  of  aerosol  transmission  (other  attenuators 
not  included)  and  show  the  excellent  agreement  obtained 
between  the  HAZE  and  LOWTRAN  results.  No  detailed  statis¬ 
tical  analysis  of  the  small  errors  is  deemed  appropriate 
because  of  the  empirical  nature  of  the  LOWTRAN  methodology. 
The  major  value  of  LOWTRAN  and  HAZE  aerosol  predictions 
lies  in  showing  in  what  regimes  aerosol  effects  are  signifi 
cant  rather  than  in  absolute  values  of  transmission.  It  is 
clear  from  the  results,  however,  that  for  even  relatively 
"clear"  atmospheres,  excluding  aerosol  effects  (as  ASDIR 
has  done)  can  introduce  errors  in  overall  transmission 
calculations  of  20  per  cent  or  more  dependent  upon  the 
aerosol  model  selected. 


LOWTRflN  4 — AEROSOL  ONLY  TRANSMISSION 
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Figure  14.  LOWTRAN/KAZE  Subroutine  Comparison  (2-6  pm) 
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Figure  15.  LOWTRAN/HAZE  Subroutine  Comparison  (2-6  pm) 
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Figure  20.  LOWTRAN/HAZE  Subroutine  Comparison  (2-6  pm) 
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Figure  21.  LOWTRAN/HAZE  Subroutine  Comparison  (2-6  pm) 
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Figure  24.  LOWTRAN/HAZE  Subroutine  Comparison  (8-12  yirt) 


Figure  25.  LOWTRAN/HAZE  Subroutine  Comparison  (8-12  ym) 


69 


Comparison  of  the  Unmodified  ASDIR  and  LOWTRAN.  To  get 
a  baseline  comparison  between  the  unmodified  ASDIR  program 
and  LOWTRAN  (without  aerosols)  several  cases  were  run  with 
the  resulting  plots  placed  on  the  same  graphs.  Figures  30- 
32  display  these  comparisons  for  two  low  altitude  cases  and 
one  high  altitude  situation.  The  "target"  in  ASDIR  was 
modeled  as  a  small  black  body  emitter  at  500  R.  j 

It  is  apparent  from  the  figures  that  ASDIR  (solid  line 
in  these  plots)  does  not  show  the  same  small  variations  in 
transmission  as  does  LOWTRAN.  This  is  because  LOWTRAN 
includes  several  absorbers  (eg.  0^,  ^0,  CH^,  N2 ,  and  0^) 
that  ASDIR  does  not  include.  ASDIR  was  written  as  a  time¬ 
wise  tool  for  preliminary  IR  engineering  design  and  thus 
neglected  those  minor  atmospheric  constituents. 

The  major  disagreement  between  ASDIR  and  LOWTRAN  lies 
in  the  region  between  A. 5  and  about  A. 8  pm.  In  this  region, 
absorption  by  nitrogen  continuum  is  significant.  LOWTRAN 
computes  a  smaller  N2  continuum  absorption  in  that  region 
than  does  ASDIR.  Resolution  of  these  differences  is  beyond 
the  scope  of  this  independent  study,  and  these  results  are 
presented  only  for  a  comparative  assessment  of  the  two 
programs . 

Addition  of  the  HAZE  Subroutine  into  ASDIR 

The  accuracy  of  the  HAZE  subroutine  in  duplicating 
LOWTRAN 's  aerosol  transmission  calculations  has  been  demon- 
srrared.  The  HAZE  subroutine  was  installed  into  ASDIR  and 


72 


COMPARISON  OP  UNMODIFIED  ASDIR  AND  LOWTRflN 
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Low  Altitude  Comparison  (Case  1)  Between  the 
Unmodified  ASDIR  and  LOWTRAN  Program  Results 


COMPARISON  OF  UNMODIFIED  RSDIR  AND  LOWTRRN 
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Figure  31.  Low  Altitude  Comparison  (Case  2)  Between  the 
Unmodified  ASDIR  and  LOWTRAN  Program  Results 
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the  results  validated  by  comparison  with  LOWTRAN.  Figure  33 
shows  the  comparisons  of  aerosol  (only)  transmission  as 
predicted  by  the  modified  ASDIR  and  LOWTRAN  programs  for  a 
short  range,  low  altitude  case  using  the  average  continental 
aerosol.  Figure  34  compares  total  transmission  plots  from 
the  two  programs.  Figures  35  and  36  are  similar  depictions 
assuming  a  maritime  aerosol  environment.  Agreement  between 
the  modified  ASDIR  and  LOWTRAN  aerosol  transmission  is 
excellent  while  the  total  transmission  shows  the  same  dif¬ 
ferences  as  in  Figs  30-32. 

Note  that  for  the  5  km  visibility,  representing  a 
"hazy"  condition,  aerosol  effects  are  seen  to  be  a  maximum 
at  the  lower  wavelengths  and  decrease  the  overall  trans¬ 
mission  by  30  per  cent  or  more  in  the  2. 0-2. 5  pm  region. 

An  analysis  of  these  types  of  results  indicate  that  on 
relatively  clear  days  (of  high  absolute  humidity)  IR  sys¬ 
tems  (eg.  FLIR)  of  equivalent  sensitivity  will  operate 
better  in  the  3-5  pm  window  while  hazy  conditions  generally 
favor  operation  in  the  8-12  pm  region.  It  is  specifically 
this  kind  of  qualitative  observation  of  system  performance 
that  the  aerosol  computations  will  help  make  when  using 
ASDIR  to  compute  target  IR  signatures. 

Concluding  Remarks 

It  is  evident  that  the  aerosol  environment  is  one  of 
the  most  important  considerations  in  the  overall  perform¬ 
ance  of  IR  systems,  especially  at  low  altitudes.  Aerosol 
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AEROSOL  ONLY  -  ASDIR (MOD. )  AND  LOWTRAN 


Modified  ASDIR  and  LOWTRAN  Comparison 
Aerosol  Transmission  (2-6  ym) 
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Figure  35.  Modified  ASDIR  and  LOWTRAN  Comparison 
Aerosol  Transmission  (2-6  ym) 
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Figure  36.  Modified  ASDIR  and  LOWTRAN  Comparison 
Total  Transmission  (2-6  pm) 


effects  must  be  taken  into  account  at  every  level  of  engi¬ 
neering  design  and  evaluation  of  thermal  imaging  devices. 
The  methods  in  this  study  to  include  aerosol  computations 
into  ASDIR  represent  a  significant  improvement  in  that 
program's  capabilities.  Additionally,  the  aerosol  treat-  ' 
ment  here  is  quite  flexible.  As  more  accurate  data  on 
aerosols  is  obtained,  this  information  can  be  incorporated 
into  ASDIR  through  simple  curve  fits  of  the  new  extinction 
coefficients.  The  resulting  analytical  expressions  can  be 
inserted  directly  into  the  BETCAL  portion  (extinction  coef¬ 
ficient  calculations)  called  by  the  HAZE  subroutine. 

Already,  the  latest  edition  of  LOWTRAN  includes  such  aerosol 
models  as  volcanic,  radiation  fog,  advection  fog,  and  high 
altitude  (meteoric)  aerosols.  At  least  some  of  these  will 
be  installed  into  ASDIR  by  this  writer  in  the  following 
months . 
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Appendix  A 

A  Short  Summary  of  Mie  Theory  and  Calculations 

The  calculation  of  scattering  and  absorption  of 
radiation  by  aerosol  particles  is  very  complex.  Several 
fundamental  assumptions  are  necessary  if  the  problem  is  to 
be  treated  at  all.  This  appendix  briefly  summarizes  the 
highlights  of  Mie  theory,  undoubtedly  the  most  widely  used 
and  experimentally  verified  tool  for  computing  extinction 
due  to  small  particles  in  the  atmosphere. 

The  original  work  by  Mie  reported  in  1908  resulted 
from  basic  research  on  scattering  by  colloidal  metal  parti¬ 
cles.  His  theory  found  early  applications  in  physical 
chemistry  but  applications  to  atmospheric  optics  were  much 
slower  to  develop.  In  its  complete  form,  Mie  theory 
describes  the  scattering  characteristics  of  a  broad  range 
of  particle  sizes  and  refractive  indices.  Additionally, 
the  theory  reduces  to  Rayleigh  theory  when  the  particles 
are  small  compared  to  the  wavelength  of  the  incident  radia¬ 
tion.  Thus,  Rayleigh  theory  is  simply  a  subset  of  the  more 
general  Mie  theory  and  adequately  describes  the  scattering 
by  the  very  small  gas  molecules  in  the  air  -  the  scattering 
which  accounts  for  the  sky's  blue  color.  Starting  with  the 
very  small  particles  (molecular  size) ,  as  the  particle  size 
relative  to  the  wavelength  increases,  there  is  a  gradual 
transition  from  Rayleigh  to  Mie  scattering.  Thus,  the 
most  important  factor  in  Mie  calculations  is  the 
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dimensionless  size  parameter  a,  defined  by 


(35) 


where  X  is  the  wavelength  in  the  medium  surrounding  the 
particle  and  r  is  the  particle  radius.  An  important  fact 
here  is  that  the  most  significant  factor,  relative  size, 
means  that  increasing  the  wavelength  always  results  in 
making  the  particle  appear  smaller.  Thus,  a  haze  particle 
with  a  radius  of  1  ym  can  be  considered  a  relatively  large 
particle  when  irradiated  by  visible  light  (approximately 
0.55  ym)  but  becomes  a  relatively  small  one  for,  say, 
energy  at  10.6  ym. 

Next,  one  should  note  that  Mie's  theory  solved  the 
problem  of  a  monochromatic  plane  wave  incident  on  a  homo¬ 
geneous,  isotropic  sphere  of  radius  r  surrounded  by  a 
transparent  homogeneous  and  isotropic  medium.  Thus,  in 
using  the  theory  we  implicitly  assume  that  the  atmospheric 
particles  are  spherical  and  isotropic.  That  assumption, 
employed  in  this  th°sis  by  the  use  of  LOWTRAN  extinction 
coefficients  calculated  by  Mie  theory,  is  quite  valid  for 
wet  haze  particles  such  as  fog  and  cloud  droplets  but 
certainly  less  valid  for  dust-like  particles  with  their 
random  shapes.  However,  scattering  by  non-spherical  parti¬ 
cles  is  extremely  complex  and  this  writer  has  encountered 
no  practical  applications  of  that  case  in  the  literature. 

An  atmospheric  particle  represents  an  optical 
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discontinuity  to  a  light  wave  because  the  particle's 
refractive  index  is  greater  than  that  of  the  surrounding 
air.  If  the  particle  only  scatters  energy,  with  no  absorp¬ 
tion  of  the  incident  flux,  then  the  index  of  refraction  is 
a  purely  real  number  from  about  1.33-1.60  (McCartney,  1976: 
224).  When  absorption  is  significant,  the  index  of  refrac¬ 
tion  m(A)  is  represented  as  a  complex  number  which  can  be 
written  as 


m(A)  =  n^(A)  -  i^(A)  (36) 

where  n^  and  ^  are  the  real  and  complex  parts,  respectively. 

The  index  of  refraction  of  particles  is  a  very  impor¬ 
tant  parameter  in  atmospheric  optics  since  it  determines, 
in  effect,  the  scattering  and  absorption  cross  sections  for 
the  aerosols.  In  the  following  section,  the  use  of  the 
size  parameter  and  the  index  of  refraction  m  to  calculate 
the  attenuation  coefficients  by  Mie  theory  will  be  outlined. 
Much  of  this  treatment  can  be  traced  back  to  Stratton's 
(1941)  or  Van  de  Hulst's  (1957)  work,  and  the  interested 
reader  is  referred  to  those  sources  for  more  detailed 
development  of  the  calculations  for  several  other  categories 
of  radiation  (eg.  polarized  energy).  McCartney  (1976)  also 
gives  an  excellent  summary  of  the  theory  and  much  of  his 
treatment  is  condensed  here. 

The  objective  of  these  calculations  is  to  compute  an 
extinction  coefficient  for  use  in  the  Beer's  law  formulation 


86 


(37) 


Tfl  ■=  exp(-BL) 

where 

=  transmittance  due  to  aerosols 
L  =  equivalent  path  length  of  radiation 

To  accomplish  this  aim,  it  is  first  noted  that  the  angular 
characteristics  of  Mie  scattering  for  all  particle  sizes 
and  wavelengths  are  expressed  by  two  intensity  distribution 
functions.  The  light  scattered  at  an  observation  angle  0 
by  a  particle  may  be  treated  as  consisting  of  two  compo¬ 
nents  having  intensities  Ij.  (0)  and  IM  (0).  These 
components  are  proportional  to  two  intensity  distribution 
functions  i^  and  ^  respectively.  These  functions  repre¬ 
sent  the  heart  of  Mie  theory  and  depend  on  the  size 
parameter  a,  the  complex  index  of  refraction  m,  and  on  the 
scattering  angle  0.  For  a  spherical,  isotropic  particle 
the  intensity  functions  may  be  expressed  as 


i,(a,m,0)  =  |  l 
n=l 


fo—  (a  it  +  b  T  ) 
n(n+l)  n  n  n  n' 


(38-a) 


l2(a,tM>  -  | <Vn  +  Vn>l‘  <38-b> 

where  the  n's  are  positive  integers;  an  and  bn  are  found 

from  Ricatti-Bessel  functions.  The  factors  a  and  b  have 

n  n 

arguments  formed  from  the  particle  characteristics  a  and  m, 


but  they  are  independent  of  the  angle  0 .  The  functions 
and  in  Eqs  38  depend  only  on  the  angle  0  and  involve  the 
first  two  derivatives  of  Legendre  polynomials  having  order 
n  and  argument  cos  0.  Extensive  tabulations  of  these  func¬ 
tions  are  formulated,  especially  since  the  advent  of  high 
speed  computers,  and  are  published  in  many  sources. 

When  a  particle  is  illuminated  by  unpolarized  light, 
the  scattered  radiation  consists  of  two  incoherent  compo¬ 
nents  according  to 


1(0)  -  E, 


fh  +  *2 


4^ : 


(39) 


where  i^,  i ^  are  defined  in  Eqs  38  and  Eq  is  the  irradiance 
of  the  incident  light.  Additionally,  the  angular  cross 
section  <jp(e)  (also  referred  to  as  a  particle  scattering 
efficiency  factor)  may  be  defined  as 


op(0) 


il  +  i2 


(40) 


This  factor  <^(0)  is  defined  as  that  area  such  that  the 
power  flowing  across  it  is  equal  to  the  scattered  power  per 
steradian  at  an  observation  angle  0.  Equation  40  again 
applies  only  to  unpolarized  incident  radiation;  different 
relations  written  in  terms  of  the  polarized  intensity 
components  hold  for  polarized  light. 

Next,  the  total  scattering  cross  section  Op  is  defined 
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as  that  cross  section  of  an  incident  wave  acted  upon  by  the 
particle  having  an  area  such  that  the  power  flowing  across 
it  is  equal  to  the  total  power  scattered  in  all  directions. 
Thus,  for  a  single  particle 

Att 

o  *  /  o  (6)  du>  (41) 

P  o  p 

where  dw  is  a  differential  increment  of  steradian.  Analo¬ 
gous  cross  sections  can  be  defined  for  polarized  incident 
light  while  Eq  41  holds  for  unpolarized  radiation.  Then, 
since 


dcj  *  2  it  sin  6  d0  (42) 

we  may  perform  the  integration  and  define  an  efficiency 
factor  Qsc,  which  is  the  ratio  of  the  scattering  cross 
section  to  the  geometric  cross  section,  as 


ir 


'•sc 


4* -4 

irr2  rz 


or  (0 )  sine  de 


(43) 


The  factor  Qgc  defines  the  "efficiency"  with  which  the 
particle  totally  scatters  the  incident  light.  From  Eq  40 
we  see  that 


Q  -  ~ 

X8C  -2 


2 

ir 

’  2lL 

\h  +  4] 

r2 

a 

4tt  2 

2 

sine  de 


(44-a) 
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(44-b) 


Q,c  -  IT  /  »1  +  *2>  sln6  de 
u  o 

where  has  previously  been  written  as  the  dimensionless 
size  ratio. 

The  integration  in  Eqs  44  is  largely  illustrative 

since  for  computational  purposes  i^  and  i ^  can  be  defined 

in  terms  of  the  complex  functions  an  and  bn  as  shown  in 

Eqs  38.  In  those  terms  a  can  be  written  as 

P 

a  =  AI  I  (2n+l)(|a  |2  +  |b  |2)  (45) 

2tt  n-l  n  n 

and  the  efficiency  factor  becomes 

Both  of  these  forms  are  well  suited  to  calculation  on 
digital  computers. 

In  a  method  completely  analogous  to  the  previous 
development  an  extinction  efficiency  factor  may  be  computed 
when  the  index  of  refraction  has  an  imaginary  part,  ie. 
when  the  particle  absorbs  energy.  The  result  is 


<W  -  b  l  (2n+D^(«n  +  V 

u  n-1 

where  Re  denotes  the  real  part  of  the  sum  of  an  and  bR. 
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Equations  46  and  47  compute  the  efficiency  factors 
for  a  single  particle.  The  extinction  due  to  N  identical 
particles  of  radius  r  in  a  volume  is  then 


6 


ext 


Nirr2  Q 


ext 


(48) 


where  the  area  is  nr2.  Generally,  as  was  discussed  in 
Chapter  II,  there  are  many  different  sized  particles  in  a 
volume,  and  thus  the  extinction  coefficient  for  a  collection 
of  polydispersed  spheres  is 

r2 

8ext  "  /*Qext(r)  n(r)  r2{ir  (49) 

rl 

where 

n(r)  *  the  size  distribution  of  the  spheres 
Qext(r)  ■  extinction  efficiency  at  radius  r 


As  shown  in  Chapter  II  the  function  n(r)  for  aerosols 
is  usually  a  smoothly  varying  one  which  may  be  described  by 
a  power-law  distribution.  Using  the  methods  outlined  here 
it  can  be  understood  how  an  aerosol  size  distribution  can 
be  assumed  and,  for  a  given  index  of  refraction,  extinction 
coefficients  can  be  computed.  This  represents  the  method¬ 
ology  employed  by  the  authors  of  the  LOWTRAN  code  and 
adapted  for  use  in  this  independent  study. 
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Appendix  B 

Derivation  of  Expressions  for  Equivalent  Path  Lengths 

To  duplicate  LOWTRAN 's  results  several  values  of 
AV(z,VIS)  were  obtained  in  successive  runs  of  the  LOWTRAN 
program.  AV(z,VIS)  is  described  in  Ref  11  as  the  equiva¬ 
lent  sea  level  path  length  for  aerosol  extinction  for 
vertical  atmospheric  paths.  Use  of  the  values  of  AV(z,VIS) 
at  each  endpoint  of  a  slant  path  is  shown  in  Eq  34  where 
they  are  used  to  calculate  the  overall  equivalent  path 
length  for  the  slant  transmission  path.  AV  values  calcu¬ 
lated  by  LOWTRAN  are  functions  of  both  visibility  and 
altitude  for  heights  below  5  km  through  the  dependence  of 
the  particle  density  upon  visibility.  Above  5  km  the 
assumed  particle  densities  are  independent  of  visibility. 

Figure  37  plots  the  results  of  LOWTRAN  calculations 
for  three  visibilities  from  0-7  km.  As  it  appears  in  the 
figure  there  is  a  linear  portion  of  each  curve  for  constant 
VIS  from  z“0  to  z«3.  In  the  region  3<zf7  the  constant  VIS 
curves  are  approximated  by  a  parabola.  All  curve  fits  in 
this  appendix  were  found  using  a  HP- 9 7  computer  curve  fit 
package  capable  of  applying  exponential,  power,  or  loga¬ 
rithmic  fits  to  sets  of  data.  First  the  z-0  intercept  was 
derived  from  a  curve  fit  of  AV(0,VIS).  The  resulting 
expression  was  found  as  a  power  law  fit: 

AV(0 , VIS)  -  25.952  VIS"1-037  (50) 
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Altitude  (km) 


Figure  37.  LOWTRAN  Values  of  Equivalent  Sea  Level 
Path  Lengths  for  Vertical  Paths  (AV(z)) 
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Next,  an  expression  for  the  slope  (SLP)  of  the  con¬ 
stant  VIS  curves  in  the  linear  region  (zS3)  was  found. 

Here  a  logarithmic  curve  fit  given  by  Eq  51  produced  a 
good  fit  to  the  data. 

SLP  =  -1.187  +  0.106  In  (VIS)  (51) 

Combining  the  two  previous  expressions  gives  a  func¬ 
tion  for  computing  AV(z,VIS)  in  the  linear  region  given  by 
Eq  52. 

AV(z, VIS)=AV(0 , VIS)  exp (SLP  x  z) 

=25.952  VIS-1' 037exp ((-1. 187+0. 1061n  VIS)z) (52) 

In  the  region  3<z<7  a  parabolic  approximation  was 
developed  from  the  standard  formula  for  a  parabola  where 
h,k  are  coordinates  of  the  parabola's  vertex: 

(y-k)  =  c(x-h) 2  (53) 

In  the  present  case  the  parabola's  vertex  is  taken  as 
(7, lnAV(7,VIS))  or  (7,0.021).  The  constant  c  is  found  for 
each  curve  (given  by  VIS)  by  matching  the  parabola  to  the 
straight  lines  at  AV(3,VIS).  The  resulting  expression  for 
AV(3,VIS)  computed  using  the  linear  region  formula  (Eq  52) 
is 

AV(3,VIS)»(25.952  VIS-1' 036) exp (-3. 559+0. 312  In  VIS)  (54) 
Substituting  these  expressions  into  the  general  equation  for 


94 


the  parabola  a  simplified  expression  for  AV(z,VIS)  in  the 
parabolic  region  is 


[ 7 ■ 0-z| 

AV(z,  VIS)»0. 0208  f— I  4  - 

0.0208 


(3<z<7) 


For  altitudes  above  7  km  the  following  functions  were 
found  to  give  satisfactory  approximations  for  AV(z) : 

7.0<z<18.5  km  AV(z)  -  10(-0.0148z  +  1.549) 

18. 5<z520. 0  km  AV(z)  =  i0(-°-136z  +  0-692) 

20.0<z  AV(z)  assumed  zero  (56) 
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P rograro  Listing  of  Subroutines  HAZE,  AVCAL,  and  BETCAL 


SUBROUTINE  HAZE(R,Z1,  72, MODEL , VIS, NU, TAUH) 

COMMON  IF  LAG 

REAL  LAM  -  -  - 

•C«*****MOOrL  =  l  FOR  LOWli'AN  fVG.  CONTI  NFNTAL  AEROSOL  KODEL  * 

C  4  MOOT  L-Z  fC k  LOHffAN  HARlTlHE  AfFCSOL  HODEL  * 

C  »NU=EREOUENCY  (C--1)  * 

ACCMf  ,n,C,0,7)  -r  .330*  f&*  it  M-B*  ?>  -C  IP.’1  (-D'Z)) 

t'CON(/.,!',C,n,7I  =  (G.  r  (A*  10.**  (-0'  7))  -23.  •  <C<  1C.4  4  (-O'Z)l  >/16. 

IF  (VIS. LI  .3.L)  MLOF  L  =  2 

L  AM=  1  . G  EO  4/NU 

IFIA0=1FL AG+1 

IF(IFLAG.NE.l)  GO  TO  SCO 

IF<VIS.GL .3.0)  GO  TO  C 

PRINT4  (•ttMtUtlHIM'  M  •  IH-  .MtimHM” 

PRIM  4 ,"* LOW  ViS)‘HLITv  INPICAT FS  POSSIBLE  FOG  CONDITIONS  *  " 
PRINT  4  , r  L  0W1  RAN  MA  RIT  j  ME  AEROSOL  MODEL  BEING  USEO  *” 

PRINT-, "  *«*•**•»>*.».«  *4*4*':  ... 

5  CONTINUE 
PNO- 2E 23 

C**'0ETEEU'1NF  IF  THE  PITH  IS  H0RI70N1AL  OR  SLANT 
IF (At  S(Zl-7?),GT.f.G)  fO  TO  ler 

C*44C0I1PUU  AH7  FOR  IMF  HORIZONTAL  PATH  . . 

IF(7l.OE.0.0.AN1.zi.LF.3.0>  GO  TO  10 
IF(7l .GT. 3.0, ANO. 71 .LE.4.0)  GO  TO  20 
IF (71 .GT. 6.0)  GO  TO  40 
lF(21.r>T.R.O.AND.7l.ur.S>0)  GO  TD  30 
10  AC=A(.(.N(  13790, 0.:.  37 , 28  78 , 0.  36  ,  Zl  > 

eC  =  OCCK(137dt,0.i  37, 2326, 0.36, 71)  -•  --  -  .  .  . . 

GO  TO  10 

20  AC  =  AC ON (137 00, O.i  37 , 1C  29  »  ,0.277 ,71) 

OC=BCON(i 3780,0. A  37 , IP 29 , 0 . 277 ,  "M ) 

GO  TO  !  0 

30  AC=ACO!.(1  378C,U.:  37,  a  27.  ,  0.  1 37 , ’l) 

-  ...  BC  =  0CON(1376O.,O.<  37, (.22, 0.137, 71)  . . 

GO  TC  60 

40  IF(71.LE.f».>  GO  TO  H 

IF (71. GT. 0.0. AND. 71. LE. 18.0)  GO  TD  42 

IF ( Z 1 . GT . 18 . C . AND .71 .L  F • 30. )  GO  TO  43 

»  AH7=C  •  0 

GO  TO  500  . 

41  FN=‘-  22. *10.**  (-.137*71) 

GO  TO  f.0  .  .... 

42  PN=7fc .t*10**(-0, 0139*71) 

GO  TO  60  --  . 

43  PN  =  3EP1.‘  10.**  (-0.137*71) 

j.  -GO  TO  GO  ...  - . . 

50  PN=AC/VIS-RC 

60  AH7=FN<  R/PNO  ...  --- 

AH7=AHZ* 0.935 

PRINT*, "THE  PATH  IS  APPROXIMATELY  HORIZONTAL" 

PRINT4, “THE  EPL  IS  ",AHZ 
GO  TO  S0C 
ICO  CONTINUE 

C**  *  COMPUTE  THE  EQUIVALENT  SEA  LEVEL  PATH  LENGTH  TOR  SLANT  PATH 
PRINT*, "THE  PATH  IS  A  SLANT  PATH" 

CALL  AVCAL (71 , VIS , AV) 

AVi=A V 

PRINT*, "71=  ",Z1,"  A Vl( KM-1) =  ",AV1 

CALL  AVCAL (Z2,VIS,AV) 


A  V2=  A  V 

PRINT', "72=  ",72,"  A  V?( KM-i >  =  ",AV2 
AH7- A  BS< ( AV1-AV?)  *R/<  ? 1  - 7  H) ) 

IF(Zl.GT.6.C.ANr).Z?.GT.D.0)  GO  TO  20T 
CFAf  =  -G.  0<*3L*  (VIS1  >  21  4  2.397*  VIS-13. 52 
CFAC-CFAC/10G. 0 
CF  AC=  1 . 0 ♦  t'F A  C 

AH7-AH7' CF AC  ... 

GO  TO  300 

200  AH7-(AVHAVD/2.G 

303  PR  I  NT  '  ,  “T  HI  ETUI  VALENT  PATH  LENGTH  IS  M,AH7 
5  CO  f  NTT  NUE 

C**"o.HPUTf  TAIJH,  THE  TRANSMISSION  AT  HAVELFNGTH  LAP 

Cf.-MuTr  THE  E  X  f !  t-'C  T  1  0!J  CO£  FFlCTEt.T  BY  A  TALL  ’0  BETCAL 
CALL  OETCAL (t  AM, MOO EL , "E  T) 

C*»  «  *■  C(iMI  UT  E  THE  OMICAL  T.'PTH 

cD-erv  ah7 

c***  *  COMI  IITF  TAUH 

1AUH=IXP(-0P> 

RE  TUB N 
ENO 


{♦«  **  >«<**♦»(•!•  '  »»*«•*  *'  * 

SUBROUTINE  AVCAL ( 7, VIS, AV) 

C  'THIS  SUBROUTINE  COMPUTES  THE  EQUIVALENT  Sfl  LEVEL  PA^H  LENGTH 
C  *  ( PEP  KM)  FOP  SLANT  PATHS.  REFER  TO  FIG.  It  IN  OPTICAL  PROPERTIES  * 

C  *  OF  '(HE  ATMOSPHERE  BY  HCCLATCHEY,  FT  AL  (13'‘2).  AV(Z)  FOR  Z  LESS  * 

C  THAN  7  KM  HAS  BEEN  DERIVED  BY  SYSTEMATIC  CJf-VE  FITS  OF  LOWTRAN 
C  'COMPUTED  VALUES  AS  A  FUNCTION  OF  BOTH  VIS  MOD  ALTITUDE  (7) 

Q  III  M*  IM.imilHIIItlll 

IF (7.LE.3.3)  GO  TO  13 
IF(7.GT.3.0.AMD.7.LE.7.0)  GO  TO  20 

TFiZ.GT.T  ,0  .AN0.7  .LE.  If  .35)  GO  TO  33  ..  .. 

IF(7.GT.18,3E.AMD.7.LC.20.C)  GO  TO  48 
IF (Z.GT  .20.0)  A  V- 0 . 0 
RETURN 

10  CONTINUE 

C»*'COMPUIt  THE  VALUE  OF  AV(f,VIS> — THE  Y  INTERCEPT  OF  AV 
A.Vl=2*  ,9617GC  <  VIS'*  (-1.03EE-84) 

C**'COMPUTF  THE  SLO°E  (SIP)  OF  T-*E  AV  CURVE  AS  rUNCTION  OF  VIS 
SLP--1 .1.CE663  ♦  0. 10SS718‘ ALGG (VIS) 

C*’  COMPUTF  AV(7,VIS)  III  THE  LINEAR  REGION 
AV=AV(.'EXP(SIP'7) 

RETUI  N 

20  CONTINUE 

C  "  ’  COMPUTE  VALUE  OF  AV(3,VTS> 

AV3= (20.661766* VI Sv  ' (-1. 03( 584) ) *  EXP (-3. E > SB 7 E ♦ . 31 2134* A LOG ( VIS > > 
C"  '  COMcUT  F  A V ( 7 , VI S I  IN  THF  PARABOLIC  F  EG  I  ON 
POV.’=<  (7. 0-7)  74.1)'“  2 
AV=0.t20f * <(AV3/. T2TE)  * (POW)) 

RETURN 

30  A  V=1 0 . 0**  (-0.01477* 7  +  1,549) 

RETUF  N 

40  AV=lt.C*' (-0.136* Z+3.fcF2) 

RETURN 

ENO 


I! 
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SUBROUTINE  BE  T  CAL  OIL  AM ,  Mj  PET) 

C  »***-•  4  4  «  ¥  4  »>  »  CtO  MC  r  I  '  V  >  *  '<*•,' 

c  *  nn cal  computes  kip:  extinction  comiciE'iT  as  a  function  of 

C  *  HAVLLLNGTII  LAM  rof<  TWO  OF  THE  LOUT  RAN  A  1  1 0S  PHrf>I  C  AEROSOL 

C  *  MODELS.  M-3  USIS  THE  AVERAGE  CONTINENTAL  MCO'L.  M-2  USES 

C  "  THE  MARITIME  MODEL.  THE  MARITIME  MODEL  SHOUL1'  NORMALLY  BE 

C  *  USED  FOR  LGH  VISIBILITY  CONDITIONS  (VIS  LtSS  I II AII  f>  KM)  SINCE 

C  *  FOG  CONDITIONS  MAY  BE  PRESENT. 

COKHON  IFLAG 
FEAL  LAM 
LAK-FLAM 

IF(H.E0.2>  GO  TO  10 

J  F  (I  r  LAG.  CO.  1 )  PRI  NT  *  ,  '*  USING  AVG.  CONIINT'ITAL  MODEL" 

IF  (LAK.GE.O.CG.ASn.LAtl.LE.2.  C)  BCT  =  -  0 . 00  3-S »  LAl  +  0  .  1C2999 
IF  ( L/  M.GT  .2.0.  AID.  LM1.LT  .7.2)  BET  =•  C.  0  012'  .  AM*  *  (2  . )  -0 . 0 1070  *  LAM 
C  4  0  «  0 1  7 DAS 

IFILfM.GE.  7. 2. AND.LAH.IT. 7.9)  R£T=-7 . 0E-'.  3’  LAM+S  .9  AGE- 02 
IF (LAM.GE.7.9.AN0.L4M.I  E.  0.  2)  PEI  =  8.  33E-0'.  M.AM  +  1.2EG7E-03 
IF (LA  M.GT  . 0.2. AND. LAM. LE. 9.0)  OFT= -S . 7 9E -0 3* L AH'  •  (2.) 

C  +  0 .17 UGE  LAM-0. Or  17$ 

IF  (LAM.GT  .9.C.ANri.LAH.LE.9.?>  B  E  T  =  0 »  023C9 

IF  (LAM.GT  .9.2.  AN!)  .LAM.  I  E.  1A  .0)  BET  =  0 . 0 1C  5  5  7  - 1 .  V-  1 E- 0  A*  L  AH 
0+136.938'  EXP(-LAH) 

I F  ( L  A  M  .  GT . 1 A . 0 )  BET=0.ri0 
RETUF.N 

•10  COK13NUE 

C*¥*  ♦♦MAKIT  1  Ml  AEP.030L  MODEL  EXTINCTION  CDEFI 1  Cl rNTS 

IF(LAH.GE.0.6&,  AN<Y.LAM.LE.?.r  )  '’FT=-  0 . 02  1  •  LAM  +  0.1E7A55 

IF (LAM. GT . 2.9. AMD . L AM. LE . 3.  v )  BET=T, 099 D 

IF  (LAM.GT  .3.9.  A  JO. LAM. I  E. 8.  2)  BET  =  0.  C  0  22) »  LAM  M2.) 

.  C-C  •  03f  1<M  LAM+3 . 2‘.  A  3G1 

IF  (LAM.GT  .3. 2.  AND.  LAM.  I  T.9.0)  BET  =  -0  .  0  IE  13  3’ LAM  -  M2.) 

C+0, 17739  LAH-0.733G59 

IF (LAM.GE. 9.0. A NO.LAM.LT. 11.0)  BET  =  0 . 0 02?  23‘ LAN"  M  2. )  - 0. 0612K  L AM  +  t 
C.3f  C<: 

IF( LAM.GE. 11.0.  AH O.LAM.LE.IA.O)  OET-C  .O&Lf  (  ‘LAM-0. 0i‘J9C> 

IF  (LAM.GT  alL.0)  f>LT  =  0  • T39 

IF  (IFLAG,  EO.l)  PR)  NT ’ t "MODEL  =  2  (MAPIT1ME  AO.jjSOI.  MODEL)" 

RETURN 

END 
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